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Description 

BIOINFORMATICALLY DETECTABLE 
GROUP OF NOVEL REGULATORY VIRAL 
AND VIRAL ASSOCIATED 
OLIGONUCLEOTIDES AND USES THEREOF 

CROSS REFERENCE TO RELATED APPLICATIONS 

[0001] This application is a continuation in part of and claims 

priority from the following patent applications, the disclo- 
sures of which applications are all hereby incorporated 
herein by reference: U.S. Patent Application Serial No. 
10/604984 filed 29-Aug-03, U.S. Patent Application Se- 
rial No. 10/604945 filed 27-Aug-03, and U.S. Patent Ap- 
plication Serial No. 10/604942 filed 27-Aug-03. This ap- 
plication also claims priority from International Applica- 
tion Number: PCT/IL 03/009998, filed 26-Nov-03, the 
disclosure of which application is hereby incorporated 
herein by reference. All of the aforesaid patent applica- 
tions are entitled "Bioinformatically Detectable Group of 



Novel Viral Regulatory Genes and Uses Thereof; This ap- 
plication also is a continuation in part of and claims prior- 
ity from the following patent applications: U.S. Patent Ap- 
plication Serial No. 10/605838 filed 30-Oct-03, and U.S. 
Patent Application Serial No. 10/604944 filed 28-Aug-03. 
Both of the aforesaid patent applications are entitled 
"Bioinformatically Detectable Group of Novel HIV Regula- 
tory Genes and Uses Thereof; This application also is a 
continuation in part of and claims priority from the fol- 
lowing patent applications: U.S. Patent Application Serial 
No. 10/605840 filed 30-Oct-03, and U.S. Patent Applica- 
tion Serial No. 10/604943 filed 28-Aug-03. Both of the 
aforesaid patent applications are entitled "Bioinformati- 
cally Detectable Group of Novel Vaccinia Regulatory Genes 
and Uses Thereof; International Application Number: 
PCT/IL 03/00998, filed 26-Nov-03, entitled "Bioinformati- 
cally Detectable Group of Novel Viral Regulatory Genes 
and Uses Thereof" claims priority from the following 
patent applications, the disclosures of which applications 
are all hereby incorporated herein by reference: U.S. 
Patent Application Serial No. 10/604984 filed 29-Aug-03, 
U.S. Patent Application Serial No. 10/604945 filed 
27-Aug-03, U.S. Patent Application Serial No. 10/604942 



filed 27-Aug-03, and U.S. Provisional Patent Application 
Serial No. 60/457788 filed 27-Mar-03. All of the afore- 
said patent applications are entitled "Bioinformatically De- 
tectable Group of Novel Viral Regulatory Genes and Uses 
Thereof"; This application also claims priority from the 
following patent applications: U.S. Patent Application Se- 
rial No. 10/605838 filed 30-Oct-03, and U.S. Patent Ap- 
plication Serial No. 10/604944 filed 28-Aug-03. Both of 
the aforesaid patent applications are entitled "Bioinformat- 
ically Detectable Group of Novel HIV Regulatory Genes and 
Uses Thereof; This application also claims priority from 
the following patent applications: U.S. Patent Application 
Serial No. 10/605840 filed 30-Oct-03, U.S. Patent Appli- 
cation Serial No. 10/604943 filed 28-Aug-03, and U.S. 
Provisional Patent Application Serial No. 60/441241 filed 
17-Jan-03. All of the aforesaid patent applications are 
entitled "Bioinformatically Detectable Group of Novel Vac- 
cinia Regulatory Genes and Uses Thereof"; U.S Patent Ap- 
plication Serial No. 10/605840, filed 30-Oct-03, entitled 
"Bioinformatically Detectable Group of Novel Vaccinia 
Regulatory Genes and Uses Thereof" is a continuation of 
and claims priority from the following patent applications, 
the disclosures of which applications are all hereby incor- 



porated herein by reference: U.S. Patent Application Serial 
No. 10/604943 filed 28-Aug-03, and U.S Provisional 
Patent Application Serial No. 60/441241 filed 17-Jan-03. 
Both of the aforesaid patent applications are entitled 
"Bioinformatically Detectable Group of Novel Vaccinia 
Regulatory Genes and Uses Thereof"; This application also 
is a continuation in part of and claims priority from the 
following patent applications: U.S. Patent Application Se- 
rial No. 10/604984 filed 29-Aug-03, U.S. Patent Applica- 
tion Serial No. 10/604945 filed 27-Aug-03, U.S. Patent 
Application Serial No. 10/604942 filed 27-Aug-03, U.S 
Provisional Patent Application Serial No. 60/457788 filed 

27- Mar-03, and U.S. Patent Application Serial No. 
10/310188 filed 5-Dec-02. All of the aforesaid patent 
applications are entitled "Bioinformatically Detectable 
Group of Novel Viral Regulatory Genes and Uses Thereof"; 
This application also is a continuation in part of and 
claims priority from the following patent applications: U.S. 
Patent Application Serial No. 10/605838 filed 30-Oct-03, 
U.S. Patent Application Serial No. 10/604944 filed 

28- Aug-03, and U.S Provisional Patent Application Serial 
No. 60/441230 filed 16-Jan-03. All of the aforesaid 
patent applications are entitled "Bioinformatically De- 



tectable Group of Novel HIV Regulatory Genes and Uses 
Thereof"; U.S Patent Application Serial No.10/605838, 
filed 30-Oct-03, entitled "Bioinformatically Detectable 
Group of Novel HIV Regulatory Genes and Uses Thereof is 
a continuation of and claims priority from the following 
patent applications, the disclosures of which applications 
are all hereby incorporated herein by reference: U.S. 
Patent Application Serial No. 10/604944 filed 28-Aug-03, 
and U.S Provisional Patent Application Serial No. 
60/441230 filed 16-Jan-03. Both of the aforesaid patent 
applications are entitled "Bioinformatically Detectable 
Group of Novel HIV Regulatory Genes and Uses Thereof; 
This application also is a continuation in part of and 
claims priority from the following patent applications: U.S. 
Patent Application Serial No. 10/604984 filed 29-Aug-03, 
U.S. Patent Application Serial No. 10/604945 filed 
27-Aug-03, U.S. Patent Application Serial No. 10/604942 
filed 27-Aug-03, U.S Provisional Patent Application Serial 
No. 60/457788 filed 27-Mar-03, and U.S. Patent Applica- 
tion Serial No. 10/310188 filed 5-Dec-02. All of the 
aforesaid patent applications are entitled "Bioinformati- 
cally Detectable Group of Novel Viral Regulatory Genes 
and Uses Thereof"; This application also is a continuation 



in part of and claims priority from the following patent 
applications: U.S. Patent Application Serial No. 10/604943 
filed 28-Aug-03, and U.S Provisional Patent Application 
Serial No. 60/441241 filed 17-Jan-03. Both of the afore- 
said patent applications are entitled "Bioinformatically De- 
tectable Group of Novel Vaccinia Regulatory Genes and 
Uses Thereof; U.S Patent Application Serial 
No.10/604984, filed 29-Aug-03, entitled "Bioinformati- 
cally Detectable Group of Novel Viral Regulatory Genes 
and Uses Thereof" is a continuation of U.S Provisional 
Patent Application Serial No. 60/457788, filed 
27-Mar-03, entitled "Bioinformatically Detectable Group 
of Novel Regulatory Genes and Uses Thereof ", the disclo- 
sure of which is hereby incorporated herein and claims 
priority therefrom; and is a continuation in part of and 
claims priority from the following patent applications, the 
disclosures of which applications are all hereby incorpo- 
rated herein by reference: U.S. Patent Application Serial 
No. U.S. Patent Application Serial No. 10/604945 filed 
27-Aug-03, U.S. Patent Application Serial No. 10/604942 
filed 27-Aug-03, U.S. Patent Application Serial No. 
10/310188 filed 5-Dec-02, and U.S. Patent Application 
Serial No. 10/303778 filed 26-Nov-02. All of the afore- 



said patent applications are entitled "Bioinformatically De- 
tectable Group of Novel Viral Regulatory Genes and Uses 
Thereof"; This application also is a continuation in part of 
and claims priority from the following patent applications: 
U.S. Patent Application Serial No. 10/604944 filed 
28-Aug-03, and U.S Provisional Patent Application Serial 
No. 60/441230 filed 16-Jan-03. Both of the aforesaid 
patent applications are entitled "Bioinformatically De- 
tectable Group of Novel HIV Regulatory Genes and Uses 
Thereof"; This application also is a continuation in part of 
and claims priority from the following patent applications: 
U.S. Patent Application Serial No. 10/604943 filed 
28-Aug-03, and U.S Provisional Patent Application Serial 
No. 60/441241 filed 17-Jan-03. Both of the aforesaid 
patent applications are entitled "Bioinformatically De- 
tectable Group of Novel Vaccinia Regulatory Genes and 
Uses Thereof; U.S Patent Application Serial 
No.10/604943 filed 28-Aug-03, entitled "Bioinformati- 
cally Detectable Group of Novel Vaccinia Regulatory Genes 
and Uses Thereof" is a continuation of U.S Provisional 
Patent Application Serial No. 60/441241, filed 17-Jan-03, 
entitled "Bioinformatically Detectable Group of Novel Vac- 
cinia Regulatory Genes and Uses Thereof ", the disclosure 



of which is hereby incorporated herein and claims priority 
therefrom; This application also is a continuation in part 
of and claims priority from the following patent applica- 
tions: U.S. Patent Application Serial No. 10/604945 filed 
27-Aug-03, U.S. Patent Application Serial No. 10/604942 
filed 27-Aug-03, U.S Provisional Patent Application Serial 
No. 60/457788 filed 27-Mar-03, U.S. Patent Application 
Serial No. 10/310188 filed 5-Dec-02, and U.S. Patent Ap- 
plication Serial No. 10/303778 filed 26-Nov-02. All of the 
aforesaid patent applications are entitled "Bioinformati- 
cally Detectable Group of Novel Viral Regulatory Genes 
and Uses Thereof"; This application also is a continuation 
in part of and claims priority from the following patent 
applications: U.S. Patent Application Serial No. 10/604944 
filed 28-Aug-03, and U.S Provisional Patent Application 
Serial No. 60/441230 filed 16-Jan-03. Both of the afore- 
said patent applications are entitled "Bioinformatically De- 
tectable Group of Novel hiv Regulatory Genes and Uses 
Thereof"; U.S Patent Application Serial No. 10/604944, 
filed 28-Aug-03, entitled "Bioinformatically Detectable 
Group of Novel HIV Regulatory Genes and Uses Thereof is 
a continuation of U.S Provisional Patent Application Serial 
No. 60/441230, filed 16-Jan-03, entitled "Bioinformati- 



cally Detectable Group of Novel HIV Regulatory Genes and 
Uses Thereof", the disclosure of which is hereby incorpo- 
rated herein and claims priority therefrom; This applica- 
tion also is a continuation in part of and claims priority 
from the following patent applications: U.S. Patent Appli- 
cation Serial No. 10/604945 filed 27-Aug-03, U.S. Patent 
Application Serial No. 10/604942 filed 27-Aug-03, U.S 
Provisional Patent Application Serial No. 60/457788 filed 
27-Mar-03, U.S. Patent Application Serial No. 10/310188 
filed 5-Dec-02, and U.S. Patent Application Serial No. 
10/303778 filed 26-Nov-02. All of the aforesaid patent 
applications are entitled "Bioinformatically Detectable 
Group of Novel Viral Regulatory Genes and Uses Thereof"; 
This application also is a continuation in part of U.S Provi- 
sional Patent Application Serial No. 60/441241, filed 
17-Jan-03, entitled "Bioinformatically Detectable Group of 
Novel Vaccinia Regulatory Genes and Uses Thereof ", the 
disclosure of which is hereby incorporated by reference 
and claims priority therefrom; U.S Patent Application Serial 
No. 10/604945, filed 27-Aug-03, entitled "Bioinformati- 
cally Detectable Group of Novel Viral Regulatory Genes 
and Uses Thereof" is a continuation of U.S Patent Applica- 
tion Serial No. 10/303778, filed 26-Nov-02, entitled 



"Bioinformatically Detectable Croup of Novel Regulatory 
Genes and Uses Thereof", the disclosure of which is 
hereby incorporated herein and claims priority therefrom; 
and is a continuation in part of and claims priority from 
the following patent applications, the disclosures of which 
applications are all hereby incorporated herein by refer- 
ence: U.S. Patent Application Serial No. 10/604942 filed 
27-Aug-03, U.S Provisional Patent Application Serial No. 
60/457788 filed 27-Mar-03, and U.S. Patent Application 
Serial No. 10/310188 filed 5-Dec-02. All of the aforesaid 
patent applications are entitled "Bioinformatically De- 
tectable Group of Novel Viral Regulatory Genes and Uses 
Thereof"; This application also is a continuation in part of 
U.S Provisional Patent Application Serial No. 60/441230 
filed 16-Jan-03, entitled "Bioinformatically Detectable 
Group of Novel HIV Regulatory Genes and Uses Thereof, 
the disclosure of which is hereby incorporated by refer- 
ence and claims priority therefrom; This application also is 
a continuation in part of U.S Provisional Patent Application 
Serial No. 60/441241, filed 17-Jan-03, entitled "Bioinfor- 
matically Detectable Group of Novel Vaccinia Regulatory 
Genes and Uses Thereof", the disclosure of which is 
hereby incorporated by reference and claims priority 



therefrom; U.S Patent Application Serial No. 10/604942, 
filed 27-Aug-03, entitled "Bioinformatically Detectable 
Group of Novel Viral Regulatory Genes and Uses Thereof" 
is a continuation of U.S Patent Application Serial No. 
10/310188, filed 5-Dec-02, entitled "Bioinformatically 
Detectable Group of Novel Regulatory Genes and Uses 
Thereof", the disclosure of which is hereby incorporated 
herein and claims priority therefrom; and is a continuation 
in part of and claims priority from the following patent 
applications, the disclosures of which applications are all 
hereby incorporated herein by reference: U.S Provisional 
Patent Application Serial No. 60/457788 filed 27-Mar-03, 
and U.S. Patent Application Serial No. 10/303778 filed 
26-Nov-02. All of the aforesaid patent applications are 
entitled "Bioinformatically Detectable Group of Novel Viral 
Regulatory Genes and Uses Thereof"; This application also 
is a continuation in part of U.S Provisional Patent Applica- 
tion Serial No. 60/441230 filed 16-Jan-03, entitled "Bioin- 
formatically Detectable Group of Novel HIV Regulatory 
Genes and Uses Thereof, the disclosure of which is 
hereby incorporated by reference and claims priority 
therefrom; This application also is a continuation in part 
of U.S Provisional Patent Application Serial No. 60/441241, 



filed 17-Jan-03, entitled "Bioinformatically Detectable 
Group of Novel Vaccinia Regulatory Genes and Uses 
Thereof", the disclosure of which is hereby incorporated 
by reference and claims priority therefrom; U.S Provisional 
Patent Application Serial No.60/457788, filed 27-Mar-03, 
entitled "Bioinformatically Detectable Group of Novel Viral 
Regulatory Genes and Uses Thereof" is a continuation in 
part of and claims priority from the following patent ap- 
plications, the disclosures of which applications are all 
hereby incorporated herein by reference: U.S Patent Appli- 
cation Serial No. 10/310188 filed 5-Dec-02, and U.S. 
Patent Application Serial No. 10/303778 filed 26-Nov-02. 
Both of the aforesaid patent applications are entitled 
"Bioinformatically Detectable Group of Novel Viral Regula- 
tory Genes and Uses Thereof; This application also is a 
continuation in part of U.S Provisional Patent Application 
Serial No. 60/441230 filed 16-Jan-03, entitled "Bioinfor- 
matically Detectable Group of Novel HIV Regulatory Genes 
and Uses Thereof", the disclosure of which is hereby in- 
corporated by reference and claims priority therefrom; 
This application also is a continuation in part of U.S Provi- 
sional Patent Application Serial No. 60/441241, filed 
17-Jan-03, entitled "Bioinformatically Detectable Group of 



Novel Vaccinia Regulatory Genes and Uses Thereof ", the 
disclosure of which is hereby incorporated by reference 
and claims priority therefrom; U.S Provisional Patent Ap- 
plication Serial No.60/441241, filed 17-Jan-03, entitled 
"Bioinformatically Detectable Group of Novel Vaccinia 
Regulatory Genes and Uses Thereof is a continuation in 
part of and claims priority from the following patent ap- 
plications, the disclosures of which applications are all 
hereby incorporated herein by reference: U.S Patent Appli- 
cation Serial No. 10/310188 filed 5-Dec-02, and U.S. 
Patent Application Serial No. 10/303778 filed 26-Nov-02. 
Both of the aforesaid patent applications are entitled 
"Bioinformatically Detectable Group of Novel Viral Regula- 
tory Genes and Uses Thereof; This application also is a 
continuation in part of U.S Provisional Patent Application 
Serial No. 60/441230 filed 16-Jan-03, entitled "Bioinfor- 
matically Detectable Group of Novel HIV Regulatory Genes 
and Uses Thereof", the disclosure of which is hereby in- 
corporated by reference and claims priority therefrom; U.S 
Provisional Patent Application Serial No. 60/441230, filed 
17-Jan-03, entitled "Bioinformatically Detectable Group of 
Novel Vaccinia Regulatory Genes and Uses Thereof is a 
continuation in part of and claims priority from the fol- 



lowing patent applications, the disclosures of which appli- 
cations are all hereby incorporated herein by reference: 
U.S Patent Application Serial No. 10/310188 filed 
5-Dec-02, and U.S. Patent Application Serial No. 
10/303778 filed 26-Nov-02. Both of the aforesaid patent 
applications are entitled "Bioinformatically Detectable 
Group of Novel Viral Regulatory Genes and Uses Thereof"; 
U.S Patent Application Serial No. 10/310188, filed 
05-Dec-02, entitled "Bioinformatically Detectable Group 
of Novel Viral Regulatory Genes and Uses Thereof" is a 
continuation in part of U.S Patent Application Serial 
No.10/303778, filed 26-Nov-02, entitled "Bioinformati- 
cally Detectable Group of Novel Viral Regulatory Genes 
and Uses Thereof ", the disclosure of which is hereby in- 
corporated by reference and claims priority therefrom. 
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Background of the invention 

FIELD OF THE INVENTION 

[0022] The present invention relates to a group of bioinformati- 
cally detectable novel viral oligonucleotides and to a 
group of bioinformatically detectable novel human 
oligonucleotides associated with viral infections, both are 
identified here as Genomic Address Messenger or CAM 
oligonucleotides. 

[0023] All of which are believed to be related to the micro RNA 

(miRNA) group of oligonucleotides. 
DESCRIPTION OF PRIOR ART 

[0024] Micro RNAs (miRNA), are short ~22nt non-coding regula- 
tory RNA oligonucleotides, found in a wide range of 
species, believed to function as specific gene translation 
repressors, sometimes involved in cell-differentiation. 

[0025] The ability to detect novel miRNAs is limited by the 



methodologies used to detect such oligonucleotides. All 
miRNAs identified so far either present a visibly discern- 
able whole body phenotype, as do Lin-4 and Let-7 
(Wightman,B., Ha,l., and Ruvkun.G., Cell 75:855-862 
(1993); Reinhart et al. Nature 403: 901-906 (2000)), or 
produce sufficient quantities of RNA so as to be detected 
by the standard molecular biological techniques. 
[0026] studies reporting miRNAs (Lau et al., Science 

294:858-862 (2001), Lagos-Quintana et al., Science 294: 
853-858 (2001)) discovered 93 miRNAs in several 
species, by sequencing a limited number of clones (300 
by Lau and 100 by Lagos-Quintana) of small segments 
(i.e. size fractionated) RNA. MiRNAs detected in these 
studies therefore, represent the more prevalent among 
the miRNA oligonucleotide family, and can not be much 
rarer than 1% of all small ~20nt-long RNA oligonu- 
cleotides. 

[0027] The aforesaid studies provide no basis for detection of 

miRNA oligonucleotides which either do not present a vi- 
sually discernable whole body phenotype, or are rare (e.g. 
rarer than 0.1% of all size fractionated ~20nt-long RNA 
segments expressed in the tissues examined), and there- 
fore do not produce significant enough quantities of RNA 



so as to be detected by standard biological techniques. To 
date, miRNAs have not been detected in viruses. 
[0028] The following U.S. Patents relate to bioinformatic detec- 
tion of genes: U.S Patent No. 348935, entitled "Statistical 
algorithms for folding and target accessibility prediction 
and design of nucleic acids", U.S Patent No. 6,369,195, 
entitled "Prostate-specific gene for diagnosis, prognosis 
and management of prostate cancer", and U.S Patent 
No. 6, 291, 666 entitled "Spike tissue-specific promoter", 

each of which is hereby incorporated by reference herein. 
BRIEF DESCRIPTION OF SEQUENCE LISTING, LARGE TABLES 
AND COMPUTER PROGRAM LISTING 

[0029] a sequence listing is attached to the present invention, 

comprising 399737 genomic sequences, is contained in a 
file named SEQ_LIST.txt (63481KB, 28-Mar-04), and is 
hereby incorporated by reference herein. 

[0030] Large tables relating to genomic sequences are attached 
to the present application, appear in 10 table files (size, 
creation date), incorporated herein: TABLEl.txt (8737KB, 
28-Mar-04); TABLE2.txt (165851KB, 28-Mar-04); TA- 
BLE3.txt (1443KB, 28-Mar-04); TABLE4.txt (8120KB, 
28-Mar-04), TABLE5.txt (1567KB, 28-Mar-04), TABLE6.txt 
(42273KB, 28-Mar-04) and TABLE7.txt (76926KB, 



28-Mar-04), TABLE8.txt (198035KB, 28-Mar-04), TA- 
BLE9.txt (30207KB, 28-Mar-04) and TABLE10.txt (5084KB, 
28-Mar-04), all of which are incorporated by reference 
herein. 

[0031] a computer program listing of a computer program con- 
structed and operative in accordance with a preferred em- 
bodiment of the present invention is enclosed on an elec- 
tronic medium in computer readable form, and is hereby 
incorporated by reference herein. The computer program 
listing is contained in 6 files, the name, sizes and creation 
date of which are as follows: AUXILARY_FILES.txt (117K, 
14-Nov-03); EDIT_DISTANCE.txt (144K,24-Nov-03); 
FIRST-K.txt (96K, 24-Nov-03); HAIRPIN_PREDICTION.txt 
(19K, 25-Mar-04); TWO_PHASED_SIDE_SELECTOR.txt (4K, 
14-Nov-03); and TWO_PHASED_PREDICTOR.txt (74K, 

14-Nov-03). 
Summary of the invention 

[0032] The present invention relates to a novel group of bioinfor- 
matically detectable viral regulatory RNA oligonucleotides, 
which repress expression of human target genes, by 
means of complementary hybridization to binding sites in 
untranslated regions of these human target genes. It is 
believed that this novel group of viral oligonucleotides 



represents a pervasive viral mechanism of attacking hosts, 
and therefore knowledge of this novel group of viral 
oligonucleotides may be useful in preventing and treating 
viral diseases. 

[0033] Additionally, the present invention relates to a novel 

group of bioinformatically detectable human regulatory 
RNA oligonucleotides, which repress expression of viral 
target genes, by means of complementary hybridization to 
binding sites in untranslated regions of these viral target 
genes. It is believed that this novel group of human 
oligonucleotides represents a pervasive novel anti-viral 
host defense mechanism, and therefore knowledge of this 
novel group of human oligonucleotides may be useful in 
preventing and treating viral diseases. 

[0034] Furthermore, the present invention relates to a novel 

group of bioinformatically detectable human regulatory 
RNA oligonucleotides, which repress expression of human 
target genes associated with viral diseases, by means of 
complementary hybridization to binding sites in untrans- 
lated regions of these human target genes. It is believed 
that this novel group of human oligonucleotides repre- 
sents a pervasive novel host response mechanism, and 
therefore knowledge of this novel group of human 



oligonucleotides may be useful in preventing and treating 
viral diseases. 

[0035] Additionally, the present invention relates to a novel 

group of bioinformatically detectable viral regulatory RNA 
oligonucleotides, which repress expression of viral target 
genes, by means of complementary hybridization to bind- 
ing sites in untranslated regions of these viral target 
genes. It is believed that this novel group of viral oligonu- 
cleotides represents a pervasive novel internal viral regu- 
lation mechanism, and therefore knowledge of this novel 
group of viral oligonucleotides may be useful in prevent- 
ing and treating viral diseases. 

[0036] | n various preferred embodiments, the present invention 
seeks to provide improved method and system for detec- 
tion and prevention of viral diseases, which are mediated 
by the abovementioned groups of novel oligonucleotides. 

[0037] Accordingly, the invention provides several substantially 
pure nucleic acids (e.g., genomic DNA, cDNA or synthetic 
DNA) each comprising a novel CAM oligonucleotide, vec- 
tors comprising the DNAs, probes comprising the DNAs, a 
method and system for selectively modulating translation 
of known target genes utilizing the vectors, and a method 
and system utilizing the GAM probes to modulate expres- 



sion of target genes. 
[0038] The present invention represents a scientific break- 
through, disclosing novel miRNA-like oligonucleotides the 
number of which is dramatically larger than previously be- 
lieved existed. Prior-art studies reporting miRNAs ((Lau et 
al., Science 294:858-862 (2001), Lagos-Quintana et al., 
Science 294: 853-858 (2001)) discovered 93 miRNAs in 
several species, including 21 in human, using conven- 
tional molecular biology methods, such as cloning and se- 
quencing. 

[0039] Molecular biology methodologies employed by these 
studies are limited in their ability to detect rare miRNA 
oligonucleotides, since these studies relied on sequencing 
of a limited number of clones (300 clones by Lau and 100 
clones by Lagos-Quintana) of small segments (i.e. size 
fractionated) of RNA. MicroRNAs detected in these studies 
therefore, represent the more prevalent among the miRNA 
oligonucleotide family, and are typically not be much rarer 
than 1% of all small ~20nt-long RNA oligonucleotides 
present in the tissue from the RNA was extracted. 

[0040] Recent studies state the number of miRNA oligonu- 
cleotides to be limited, and describe the limited sensitivity 
of available methods for detection of miRNA: The estimate 



of 255 human miRNA oligonucleotides is an upper bound 
implying that no more than 40 miRNA oligonucleotides 
remain to be identified in mammals (Lim et al., Science, 
299:1540 (2003)); Estimates place the total number of 
vertebrate miRNA genes at about 200-250 (Ambros et al. 
Curr. Biol. 13:807-818 (2003)); and Confirmation of very 
low abundance miRNAs awaits the application of detection 
methods more sensitive than Northern blots (Ambros et 
al. Curr. Biol. 13:807-818 (2003)). 
[0041] The oligonucleotides of the present invention represent a 
revolutionary new dimension of genomics and of biology: 
a dimension comprising a huge number of non-protein 
coding oligonucleotides which modulate expression of 
thousands of proteins and are associated with numerous 
major diseases. This new dimension disclosed by the 
present invention dismantles a central dogma that has 
dominated life-sciences during the past 50 years, a 
dogma which has emphasized the importance of protein 
coding regions of the genome, holding non-protein cod- 
ing regions to be of little consequence, often dubbing 
them junk DNA. 

[0042] indeed, only in November, 2003 has this long held belief 
as to the low importance of non-protein coding regions 



been vocally challenged. As an example, an article titled 
The Unseen Genome - Gems in the Junk (Gibbs, W.W. Sci. 
Am. 289:46-53 (2003)) asserts that the failure to recog- 
nize the importance of non-protein- coding regions may 
well go down as one of the biggest mistakes in the history 
of molecular biology. Gibbs further asserts that what was 
dammned as junk because it was not understood, may in 
fact turn out to be the very basis of human complexity. 
The present invention provides a dramatic leap in under- 
standing specific important roles of non-protein coding 
regions. 

[0043] An additional scientific breakthrough of the present in- 
vention is a novel conceptual model disclosed by the 
present invention, which conceptual model is preferably 
used to encode in a genome the determination of cell- 
differentiation, utilizing oligonucleotides and polynu- 
cleotides of the present invention. 

[0044] | n various preferred embodiments, the present invention 
seeks to provide an improved method and system for 
specific modulation of the expression of specific target 
genes involved in significant human diseases. It also pro- 
vides an improved method and system for detection of the 
expression of novel oligonucleotides of the present inven- 



tion, which modulate these target genes. In many cases 
the target genes may be known and fully characterized, 
however in alternative embodiments of the present inven- 
tion, unknown or less well characterized genes may be 
targeted. 

[0045] a Nucleic acid is defined as a ribonucleic acid (RNA) 

molecule, or a deoxyribonucleic acid (DNA) molecule, or 
complementary deoxyribonucleic acid (cDNA), comprising 
either naturally occurring nucleotides or non-naturally oc- 
curring nucleotides. 

[0046] Substantially pure nucleic acid, Isolated Nucleic Acid, Iso- 
lated Oligoucleotide and Isolated Polynucleotide are de- 
fined as a nucleic acid that is free of the genome of the 
organism from which the nucleic acid is derived, and in- 
clude, for example, a recombinant nucleic acid which is 
incorporated into a vector, into an autonomously replicat- 
ing plasmid or virus, or into the genomic nucleic acid of a 
prokaryote or eukaryote at a site other than its natural 
site; or which exists as a separate molecule (e.g., a cDNA 
or a genomic or cDNA fragment produced by PCR or re- 
striction endonuclease digestion) independent of other 
nucleic acids. 

[0047] A n oligonucleotide is defined as a nucleic acid comprising 



2-139 nucleotides, or preferably 16-120 nucleotides. A 
Polynucleotide is defined as a nucleic acid comprising 
140-5000 nucleotides, or preferably 140-1000 nu- 
cleotides. 

[0048] a Complementary sequence is defined as a first nu- 
cleotide sequence which reverses complementary of a 
second nucleotide sequence: the first nucleotide sequence 
is reversed relative to a second nucleotide sequence, and 
wherein each nucleotide in the first nucleotide sequence is 
complementary to a corresponding nucleotide in the sec- 
ond nucleotide sequence (e.g. ATGGC is the complemen- 
tary sequence of CCCAT). 

[0049] Hybridization, Binding and Annealing are defined as hy- 
bridization, under in-vivo physiologic conditions, of a first 
nucleic acid to a second nucleic acid, which second nu- 
cleic acid is at least partially complementary to the first 
nucleic acid. 

[0050] a Hairpin Sturucture is defined as an oligonucleotide hav- 
ing a nucleotide sequence that is 50-140 nucleotides in 
length, the first half of which nucleotide sequence is at 
least partially complementary to the second part thereof, 
thereby causing the nucleic acid to fold onto itself, form- 
ing a secondary hairpin structure. 



[0051] a Hairpin Shaped Precursor is defined as a Hairpin Struc- 
ture which is processed by a Dicer enzyme complex, 
yielding an oligonucleotide which is about 19 to about 24 
nucleotides in length. 

[0052] "inhibiting translation" is defined as the ability to prevent 
synthesis of a specific protein encoded by a respective 
gene by means of inhibiting the translation of the mRNA 
of this gene. For example, inhibiting translation may in- 
clude the following steps: (1) a DNA segment encodes an 
RNA, the first half of whose sequence is partially comple- 
mentary to the second half thereof; (2) the precursor folds 
onto itself forming a hairpin-shaped precursor; (3) a Dicer 
enzyme complex cuts the hairpin shaped precursor yield- 
ing an oligonucleotide that is approximately 22nt in 
length; (4) the oligonucleotide binds complementarily to 
at least one binding site, having a nucleotide sequence 
that is at least partially complentary to the oligonu- 
cleotide, which binding site is located in the mRNA of a 
target gene, preferably in the untranslated region (UTR) of 
a target gene, such that the binding inhibts translation of 
the target protein. 

[0053] a "Translation inhibitor site" is defined as the minimal nu- 
cleotide sequence sufficient to inhibit translation. 



[0054] There is thus provided in accordance with a preferred em- 
bodiment of the present invention A bioinformatically de- 
tectable isolated oligonucleotide which is endogenously 
processed from a hairpin-shaped precursor, and anneals 
to a portion of a mRNA transcript of a target gene, 
wherein binding of the oligonucleotide to the mRNA tran- 
script represses expression of the target gene, and 
wherein the oligonucleotide has at least 80% sequence 
identity with a nucleotide sequence selected from the 
group consisting of SEQ ID NOs: 1-14456. A bioinformat- 
ically detectable isolated oligonucleotide having a nu- 
cleotide sequence selected from the group consisting of 
SEQ ID NOs: 1-14456. A bioinformatically detectable first 
oligonucleotide which is a portion of a mRNA transcript of 
a target gene, and anneals to a second oligonucleotide 
that is endogenously processed from a hairpin precursor, 
wherein binding of the first oligonucleotide to the second 
oligonucleotide represses expression of the target gene, 
and wherein nucleotide sequence of the second nucleotide 
is selected from the group consisting of SEQ ID NOs: 
1-14456. A bioinformatically detectable oligonucleotide 
having a nucleotide sequence selected from the group 
consisting of SEQ ID NOs: 1-14456. A bioinformatically 



detectable isolated oligonucleotide which anneals to a 
portion of a mRNA transcript of a target gene associated 
with Human immunodeficiency virus 1 (HIV-1) infection, 
wherein binding of the oligonucleotide to the mRNA tran- 
script represses expression of the target gene, and 
wherein the oligonucleotide has at least 80% sequence 
identity with a nucleotide sequence selected from the 
group consisting of: (a) a sequence selected from the 
group consisting of SEQ ID NOs: 1, 3, 5 and 
63804-66308, and (b) the complement of a sequence se- 
lected from the group consisting of SEQ ID NOs: 196685, 
197420, 198047 and 194813-200687. A bioinformati- 
cally detectable first oligonucleotide which is a portion of 
a mRNA transcript of a target gene associated with Human 
immunodeficiency virus 1 (HIV-1), and anneals to a sec- 
ond oligonucleotide that is endogenously processed from 
a hairpin precursor, wherein binding of the first oligonu- 
cleotide to the second oligonucleotide represses expres- 
sion of the target gene, and wherein nucleotide sequence 
of the second nucleotide is selected from the group con- 
sisting of SEQ ID NOs: 1, 3, 5 and 63804-66308. A bioin- 
formatically detectable isolated oligonucleotide which an- 
neals to a portion of a mRNA transcript of a target gene 



associated with Human immunodeficiency virus 2 (HIV-2) 
infection, wherein binding of the oligonucleotide to the 
mRNA transcript represses expression of the target gene, 
and wherein the oligonucleotide has at least 80% se- 
quence identity with a nucleotide sequence selected from 
the group consisting of: (a) a sequence selected from the 
group consisting of SEQ ID NOs: 1, 3, 5 and 
66309-68740, and (b) the complement of a sequence se- 
lected from the group consisting of SEQ ID NOs: 202601, 
204421, 205326 and 200688-206410. A bioinformati- 
cally detectable first oligonucleotide which is a portion of 
a mRNA transcript of a target gene associated with Human 
immunodeficiency virus 2 (HIV-2), and anneals to a sec- 
ond oligonucleotide that is endogenously processed from 
a hairpin precursor, wherein binding of the first oligonu- 
cleotide to the second oligonucleotide represses expres- 
sion of the target gene, and wherein nucleotide sequence 
of the second nucleotide is selected from the group con- 
sisting of SEQ ID NOs: 1, 3, 5 and 66309-68740. A bioin- 
formatically detectable isolated oligonucleotide which an- 
neals to a portion of a mRNA transcript of a target gene 
associated with Human adenovirus A infection, wherein 
binding of the oligonucleotide to the mRNA transcript re- 



presses expression of the target gene, and wherein the 
oligonucleotide has at least 80% sequence identity with a 
nucleotide sequence selected from the group consisting 
of: (a) a sequence selected from the group consisting of 
SEQ ID NOs: 3 and 18799-20253, and (b) the complement 
of a sequence selected from the group consisting of SEQ 
ID NOs: 99724 and 99658-102371. A bioinformatically 
detectable first oligonucleotide which is a portion of a 
mRNA transcript of a target gene associated with Human 
adenovirus A, and anneals to a second oligonucleotide 
that is endogenously processed from a hairpin precursor, 
wherein binding of the first oligonucleotide to the second 
oligonucleotide represses expression of the target gene, 
and wherein nucleotide sequence of the second nucleotide 
is selected from the group consisting of SEQ ID NOs: 3 
and 18799-20253. A bioinformatically detectable isolated 
oligonucleotide which anneals to a portion of a mRNA 
transcript of a target gene associated with Human her- 
pesvirus 1 infection, wherein binding of the oligonu- 
cleotide to the mRNA transcript represses expression of 
the target gene, and wherein the oligonucleotide has at 
least 80% sequence identity with a nucleotide sequence 
selected from the group consisting of: (a) a sequence se- 



lected from the group consisting of SEQ ID NOs: 1, 3, 4 
and 29296-31435, and (b) the complement of a sequence 
selected from the group consisting of SEQ ID NOs: 
119609, 119788 and 119381-124109. A bioinformati- 
cally detectable first oligonucleotide which is a portion of 
a mRNA transcript of a target gene associated with Human 
herpesvirus 1, and anneals to a second oligonucleotide 
that is endogenously processed from a hairpin precursor, 
wherein binding of the first oligonucleotide to the second 
oligonucleotide represses expression of the target gene, 
and wherein nucleotide sequence of the second nucleotide 
is selected from the group consisting of SEQ ID NOs: 1, 3, 
4 and 29296-31435. A bioinformatically detectable iso- 
lated oligonucleotide which anneals to a portion of a 
mRNA transcript of a target gene associated with Human 
herpesvirus 4 (Epstein-Barr virus) infection, wherein bind- 
ing of the oligonucleotide to the mRNA transcript re- 
presses expression of the target gene, and wherein the 
oligonucleotide has at least 80% sequence identity with a 
nucleotide sequence selected from the group consisting 
of: (a) a sequence selected from the group consisting of 
SEQ ID NOs: 1, 3, 4 and 36686-40656, and (b) the com- 
plement of a sequence selected from the group consisting 



ofSEQIDNOs: 134996, 135760 and 134776-142920. A 
bioinformatically detectable first oligonucleotide which is 
a portion of a mRNA transcript of a target gene associated 
with Human herpesvirus 4 (Epstein-Barr virus), and an- 
neals to a second oligonucleotide that is endogenously 
processed from a hairpin precursor, wherein binding of 
the first oligonucleotide to the second oligonucleotide re- 
presses expression of the target gene, and wherein nu- 
cleotide sequence of the second nucleotide is selected 
from the group consisting of SEQ ID NOs: 1, 3, 4 and 
36686-40656. A bioinformatically detectable isolated 
oligonucleotide which anneals to a portion of a mRNA 
transcript of a target gene associated with Human aden- 
ovirus B (HAdV-B) infection, wherein binding of the 
oligonucleotide to the mRNA transcript represses expres- 
sion of the target gene, and wherein the oligonucleotide 
has at least 80% sequence identity with a nucleotide se- 
quence selected from the group consisting of: (a) a se- 
quence selected from the group consisting of SEQ ID NOs: 
3 and 20254-22197, and (b) the complement of a se- 
quence selected from the group consisting of SEQ ID NOs: 
102401 and 102372-106136. A bioinformatically de- 
tectable first oligonucleotide which is a portion of a mRNA 



transcript of a target gene associated with Human aden- 
ovirus B (HAdV-B), and anneals to a second oligonu- 
cleotide that is endogenously processed from a hairpin 
precursor, wherein binding of the first oligonucleotide to 
the second oligonucleotide represses expression of the 
target gene, and wherein nucleotide sequence of the sec- 
ond nucleotide is selected from the group consisting of 
SEQ ID NOs: 3 and 20254-22197. A bioinformatically de- 
tectable isolated oligonucleotide which anneals to a por- 
tion of a mRNA transcript of a target gene associated with 
Human adenovirus C infection, wherein binding of the 
oligonucleotide to the mRNA transcript represses expres- 
sion of the target gene, and wherein the oligonucleotide 
has at least 80% sequence identity with a nucleotide se- 
quence selected from the group consisting of: (a) a se- 
quence selected from the group consisting of SEQ ID NOs: 
3 and 22198-23499, and (b) the complement of a se- 
quence selected from the group consisting of SEQ ID NOs: 
106234 and 106137-108504. A bioinformatically de- 
tectable first oligonucleotide which is a portion of a mRNA 
transcript of a target gene associated with Human aden- 
ovirus C, and anneals to a second oligonucleotide that is 
endogenously processed from a hairpin precursor, 



wherein binding of the first oligonucleotide to the second 
oligonucleotide represses expression of the target gene, 
and wherein nucleotide sequence of the second nucleotide 
is selected from the group consisting of SEQ ID NOs: 3 
and 22198-23499. A bioinformatically detectable isolated 
oligonucleotide which anneals to a portion of a mRNA 
transcript of a target gene associated with Human aden- 
ovirus D infection, wherein binding of the oligonucleotide 
to the mRNA transcript represses expression of the target 
gene, and wherein the oligonucleotide has at least 80% 
sequence identity with a nucleotide sequence selected 
from the group consisting of: (a) a sequence selected from 
the group consisting of SEQ ID NOs: 3 and 23500-24998, 
and (b) the complement of a sequence selected from the 
group consisting of SEQ ID NOs: 108595 and 
108505-111292. A bioinformatically detectable first 
oligonucleotide which is a portion of a mRNA transcript of 
a target gene associated with Human adenovirus D, and 
anneals to a second oligonucleotide that is endogenously 
processed from a hairpin precursor, wherein binding of 
the first oligonucleotide to the second oligonucleotide re- 
presses expression of the target gene, and wherein nu- 
cleotide sequence of the second nucleotide is selected 



from the group consisting of SEQ ID NOs: 3 and 
23500-24998. A bioinformatically detectable isolated 
oligonucleotide which anneals to a portion of a mRNA 
transcript of a target gene associated with Human aden- 
ovirus E infection, wherein binding of the oligonucleotide 
to the mRNA transcript represses expression of the target 
gene, and wherein the oligonucleotide has at least 80% 
sequence identity with a nucleotide sequence selected 
from the group consisting of: (a) a sequence selected from 
the group consisting of SEQ ID NOs: 3 and 24999-26195, 
and (b) the complement of a sequence selected from the 
group consisting of SEQ ID NOs: 111356 and 
111293-113584. A bioinformatically detectable first 
oligonucleotide which is a portion of a mRNA transcript of 
a target gene associated with Human adenovirus E, and 
anneals to a second oligonucleotide that is endogenously 
processed from a hairpin precursor, wherein binding of 
the first oligonucleotide to the second oligonucleotide re- 
presses expression of the target gene, and wherein nu- 
cleotide sequence of the second nucleotide is selected 
from the group consisting of SEQ ID NOs: 3 and 
24999-26195. A bioinformatically detectable isolated 
oligonucleotide which anneals to a portion of a mRNA 



transcript of a target gene associated with Human aden- 
ovirus F infection, wherein binding of the oligonucleotide 
to the mRNA transcript represses expression of the target 
gene, and wherein the oligonucleotide has at least 80% 
sequence identity with a nucleotide sequence selected 
from the group consisting of: (a) a sequence selected from 
the group consisting of SEQ ID NOs: 3 and 26196-27681, 
and (b) the complement of a sequence selected from the 
group consisting of SEQ ID NOs: 113596 and 
113585-116354. A bioinformatically detectable first 
oligonucleotide which is a portion of a mRNA transcript of 
a target gene associated with Human adenovirus F, and 
anneals to a second oligonucleotide that is endogenously 
processed from a hairpin precursor, wherein binding of 
the first oligonucleotide to the second oligonucleotide re- 
presses expression of the target gene, and wherein nu- 
cleotide sequence of the second nucleotide is selected 
from the group consisting of SEQ ID NOs: 3 and 
26196-27681. A bioinformatically detectable isolated 
oligonucleotide which anneals to a portion of a mRNA 
transcript of a target gene associated with Human her- 
pesvirus 2 infection, wherein binding of the oligonu- 
cleotide to the mRNA transcript represses expression of 



the target gene, and wherein the oligonucleotide has at 
least 80% sequence identity with a nucleotide sequence 
selected from the group consisting of: (a) a sequence se- 
lected from the group consisting of SEQ ID NOs: 1, 3, 4 
and 31436-33484, and (b) the complement of a sequence 
selected from the group consisting of SEQ ID NOs: 
125110, 125697 and 124110-128336. A bioinformati- 
cally detectable first oligonucleotide which is a portion of 
a mRNA transcript of a target gene associated with Human 
herpesvirus 2, and anneals to a second oligonucleotide 
that is endogenously processed from a hairpin precursor, 
wherein binding of the first oligonucleotide to the second 
oligonucleotide represses expression of the target gene, 
and wherein nucleotide sequence of the second nucleotide 
is selected from the group consisting of SEQ ID NOs: 1, 3, 
4 and 31436-33484. A bioinformatically detectable iso- 
lated oligonucleotide which anneals to a portion of a 
mRNA transcript of a target gene associated with Human 
herpesvirus 3 infection, wherein binding of the oligonu- 
cleotide to the mRNA transcript represses expression of 
the target gene, and wherein the oligonucleotide has at 
least 80% sequence identity with a nucleotide sequence 
selected from the group consisting of: (a) a sequence se- 



lected from the group consisting of SEQ ID NOs: 1, 3, 4 
and 33485-36685, and (b) the complement of a sequence 
selected from the group consisting of SEQ ID NOs: 
128645, 128879 and 128337-134775. A bioinformati- 
cally detectable first oligonucleotide which is a portion of 
a mRNA transcript of a target gene associated with Human 
herpesvirus 3, and anneals to a second oligonucleotide 
that is endogenously processed from a hairpin precursor, 
wherein binding of the first oligonucleotide to the second 
oligonucleotide represses expression of the target gene, 
and wherein nucleotide sequence of the second nucleotide 
is selected from the group consisting of SEQ ID NOs: 1, 3, 
4 and 33485-36685. A bioinformatically detectable iso- 
lated oligonucleotide which anneals to a portion of a 
mRNA transcript of a target gene associated with Human 
herpesvirus 5 infection, wherein binding of the oligonu- 
cleotide to the mRNA transcript represses expression of 
the target gene, and wherein the oligonucleotide has at 
least 80% sequence identity with a nucleotide sequence 
selected from the group consisting of: (a) a sequence se- 
lected from the group consisting of SEQ ID NOs: 1, 3, 4 
and 40657-45367, and (b) the complement of a sequence 
selected from the group consisting of SEQ ID NOs: 



143584, 143602 and 142921-153787. A bioinformati- 
cally detectable first oligonucleotide which is a portion of 
a mRNA transcript of a target gene associated with Human 
herpesvirus 5, and anneals to a second oligonucleotide 
that is endogenously processed from a hairpin precursor, 
wherein binding of the first oligonucleotide to the second 
oligonucleotide represses expression of the target gene, 
and wherein nucleotide sequence of the second nucleotide 
is selected from the group consisting of SEQ ID NOs: 1, 3, 
4 and 40657-45367. A bioinformatically detectable iso- 
lated oligonucleotide which anneals to a portion of a 
mRNA transcript of a target gene associated with Human 
herpesvirus 6 infection, wherein binding of the oligonu- 
cleotide to the mRNA transcript represses expression of 
the target gene, and wherein the oligonucleotide has at 
least 80% sequence identity with a nucleotide sequence 
selected from the group consisting of: (a) a sequence se- 
lected from the group consisting of SEQ ID NOs: 1, 3, 4 
and 45368-49938, and (b) the complement of a sequence 
selected from the group consisting of SEQ ID NOs: 
154384, 155376 and 153788-163278. A bioinformati- 
cally detectable first oligonucleotide which is a portion of 
a mRNA transcript of a target gene associated with Human 



herpesvirus 6, and anneals to a second oligonucleotide 
that is endogenously processed from a hairpin precursor, 
wherein binding of the first oligonucleotide to the second 
oligonucleotide represses expression of the target gene, 
and wherein nucleotide sequence of the second nucleotide 
is selected from the group consisting of SEQ ID NOs: 1, 3, 
4 and 45368-49938. A bioinformatically detectable iso- 
lated oligonucleotide which anneals to a portion of a 
mRNA transcript of a target gene associated with Human 
herpesvirus 6B infection, wherein binding of the oligonu- 
cleotide to the mRNA transcript represses expression of 
the target gene, and wherein the oligonucleotide has at 
least 80% sequence identity with a nucleotide sequence 
selected from the group consisting of: (a) a sequence se- 
lected from the group consisting of SEQ ID NOs: 1, 3, 4 
and 49939-54807, and (b) the complement of a sequence 
selected from the group consisting of SEQ ID NOs: 
163594, 163624 and 163279-174169. A bioinformati- 
cally detectable first oligonucleotide which is a portion of 
a mRNA transcript of a target gene associated with Human 
herpesvirus 6B, and anneals to a second oligonucleotide 
that is endogenously processed from a hairpin precursor, 
wherein binding of the first oligonucleotide to the second 



oligonucleotide represses expression of the target gene, 
and wherein nucleotide sequence of the second nucleotide 
is selected from the group consisting of SEQ ID NOs: 1, 3, 
4 and 49939-54807. A bioinformatically detectable iso- 
lated oligonucleotide which anneals to a portion of a 
mRNA transcript of a target gene associated with Human 
herpesvirus 7 infection, wherein binding of the oligonu- 
cleotide to the mRNA transcript represses expression of 
the target gene, and wherein the oligonucleotide has at 
least 80% sequence identity with a nucleotide sequence 
selected from the group consisting of: (a) a sequence se- 
lected from the group consisting of SEQ ID NOs: 1, 3, 4 
and 54808-59933, and (b) the complement of a sequence 
selected from the group consisting of SEQ ID NOs: 
175525, 175696 and 174170-187099. A bioinformati- 
cally detectable first oligonucleotide which is a portion of 
a mRNA transcript of a target gene associated with Human 
herpesvirus 7, and anneals to a second oligonucleotide 
that is endogenously processed from a hairpin precursor, 
wherein binding of the first oligonucleotide to the second 
oligonucleotide represses expression of the target gene, 
and wherein nucleotide sequence of the second nucleotide 
is selected from the group consisting of SEQ ID NOs: 1, 3, 



4 and 54808-59933. A bioinformatically detectable iso- 
lated oligonucleotide which anneals to a portion of a 
mRNA transcript of a target gene associated with Japanese 
encephalitis virus infection, wherein binding of the 
oligonucleotide to the mRNA transcript represses expres- 
sion of the target gene, and wherein the oligonucleotide 
has at least 80% sequence identity with a nucleotide se- 
quence selected from the group consisting of: (a) a se- 
quence selected from the group consisting of SEQ ID NOs: 
1, 2, 3, 4 and 76896-78717, and (b) the complement of a 
sequence selected from the group consisting of SEQ ID 
NOs: 222664, 224180, 224385, 224863 and 
221707-225968. A bioinformatically detectable first 
oligonucleotide which is a portion of a mRNA transcript of 
a target gene associated with Japanese encephalitis virus, 
and anneals to a second oligonucleotide that is endoge- 
nously processed from a hairpin precursor, wherein bind- 
ing of the first oligonucleotide to the second oligonu- 
cleotide represses expression of the target gene, and 
wherein nucleotide sequence of the second nucleotide is 
selected from the group consisting of SEQ ID NOs: 1, 2, 3, 
4 and 76896-78717. A bioinformatically detectable iso- 
lated oligonucleotide which anneals to a portion of a 



mRNA transcript of a target gene associated with Murray 
Valley encephalitis virus infection, wherein binding of the 
oligonucleotide to the mRNA transcript represses expres- 
sion of the target gene, and wherein the oligonucleotide 
has at least 80% sequence identity with a nucleotide se- 
quence selected from the group consisting of: (a) a se- 
quence selected from the group consisting of SEQ ID NOs: 
1, 2, 3, 4 and 83567-85368, and (b) the complement of a 
sequence selected from the group consisting of SEQ ID 
NOs: 240498, 240554, 241922, 243206 and 
240222-244441. A bioinformatically detectable first 
oligonucleotide which is a portion of a mRNA transcript of 
a target gene associated with Murray Valley encephalitis 
virus, and anneals to a second oligonucleotide that is en- 
dogenously processed from a hairpin precursor, wherein 
binding of the first oligonucleotide to the second oligonu- 
cleotide represses expression of the target gene, and 
wherein nucleotide sequence of the second nucleotide is 
selected from the group consisting of SEQ ID NOs: 1, 2, 3, 
4 and 83567-85368. A bioinformatically detectable iso- 
lated oligonucleotide which anneals to a portion of a 
mRNA transcript of a target gene associated with West 
Nile virus infection, wherein binding of the oligonu- 



cleotide to the mRNA transcript represses expression of 
the target gene, and wherein the oligonucleotide has at 
least 80% sequence identity with a nucleotide sequence 
selected from the group consisting of: (a) a sequence se- 
lected from the group consisting of SEQ ID NOs: 5 and 
89176-89792, and (b) the complement of a sequence se- 
lected from the group consisting of SEQ ID NOs: 251834 
and 251558-253039. A bioinformatically detectable first 
oligonucleotide which is a portion of a mRNA transcript of 
a target gene associated with West Nile virus, and anneals 
to a second oligonucleotide that is endogenously pro- 
cessed from a hairpin precursor, wherein binding of the 
first oligonucleotide to the second oligonucleotide re- 
presses expression of the target gene, and wherein nu- 
cleotide sequence of the second nucleotide is selected 
from the group consisting of SEQ ID NOs: 5 and 
89176-89792. A bioinformatically detectable isolated 
oligonucleotide which anneals to a portion of a mRNA 
transcript of a target gene associated with Western equine 
encephalomyelitis virus infection, wherein binding of the 
oligonucleotide to the mRNA transcript represses expres- 
sion of the target gene, and wherein the oligonucleotide 
has at least 80% sequence identity with a nucleotide se- 



quence selected from the group consisting of: (a) a se- 
quence selected from the group consisting of SEQ ID NOs: 
5 and 89793-90322, and (b) the complement of a se- 
quence selected from the group consisting of SEQ ID NOs: 
253274 and 253040-254154. A bioinformatically de- 
tectable first oligonucleotide which is a portion of a mRNA 
transcript of a target gene associated with Western equine 
encephalomyelitis virus, and anneals to a second oligonu- 
cleotide that is endogenously processed from a hairpin 
precursor, wherein binding of the first oligonucleotide to 
the second oligonucleotide represses expression of the 
target gene, and wherein nucleotide sequence of the sec- 
ond nucleotide is selected from the group consisting of 
SEQ ID NOs: 5 and 89793-90322. A bioinformatically de- 
tectable isolated oligonucleotide which anneals to a por- 
tion of a mRNA transcript of a target gene associated with 
Eastern equine encephalitis virus infection, wherein bind- 
ing of the oligonucleotide to the mRNA transcript re- 
presses expression of the target gene, and wherein the 
oligonucleotide has at least 80% sequence identity with a 
nucleotide sequence selected from the group consisting 
of: (a) a sequence selected from the group consisting of 
SEQ ID NOs: 1, 2, 3, 4 and 16459-18368, and (b) the 



complement of a sequence selected from the group con- 
sisting of SEQ ID NOs: 95931, 96664, 98184, 98729 and 
94518-98947. A bioinformatically detectable first 
oligonucleotide which is a portion of a mRNA transcript of 
a target gene associated with Eastern equine encephalitis 
virus, and anneals to a second oligonucleotide that is en- 
dogenously processed from a hairpin precursor, wherein 
binding of the first oligonucleotide to the second oligonu- 
cleotide represses expression of the target gene, and 
wherein nucleotide sequence of the second nucleotide is 
selected from the group consisting of SEQ ID NOs: 1, 2, 3, 
4 and 16459-18368. A method for bioinformatic detec- 
tion of microRNA oligonucleotides, the method including: 
bioinformatically detecting a hairpin shaped precursor 
oligonucleotide, bioinformatically detecting an oligonu- 
cleotide which is endogenously processed from the hair- 
pin shaped precursor oligonucleotide, and bioinformati- 
cally detecting a target gene of the oligonucleotide 
wherein the oligonucleotide anneals to at least one por- 
tion of a mRNA transcript of the target gene, and wherein 
the binding represses expression of the target gene, and 

the target gene is associated with a disease. 
Brief Description of Drawings 



[0055] pig.l is a simplified diagram illustrating a mode by which 
an oligonucleotide of a novel group of oligonucleotides of 
the present invention modulates expression of known tar- 
get genes; 

[0056] pig. 2 is a simplified block diagram illustrating a bioinfor- 
matic oligonucleotide detection system capable of detect- 
ing oligonucleotides of the novel group of oligonu- 
cleotides of the present invention, which system is con- 
structed and operative in accordance with a preferred em- 
bodiment of the present invention; 

[0057] pig. 3 is a simplified flowchart illustrating operation of a 
mechanism for training of a computer system to recog- 
nize the novel oligonucleotides of the present invention, 
which mechanism is constructed and operative in accor- 
dance with a preferred embodiment of the present inven- 
tion; 

[0058] Fig. 4A is a simplified block diagram of a non-coding ge- 
nomic sequence detector constructed and operative in ac- 
cordance with a preferred embodiment of the present in- 
vention; 

[0059] Fig. 4B is a simplified flowchart illustrating operation of a 
non-coding genomic sequence detector constructed and 
operative in accordance with a preferred embodiment of 



the present invention; 

[0060] pig. 5A is a simplified block diagram of a hairpin detector 
constructed and operative in accordance with a preferred 
embodiment of the present invention; 

[0061] pig. 5B is a simplified flowchart illustrating operation of a 
hairpin detector constructed and operative in accordance 
with a preferred embodiment of the present invention; 

[0062] Fig. 6A is a simplified block diagram of a dicer-cut loca- 
tion detector constructed and operative in accordance 
with a preferred embodiment of the present invention; 

[0063] Fig. 6B is a simplified flowchart illustrating training of a 
dicer-cut location detector constructed and operative in 
accordance with a preferred embodiment of the present 
invention; 

[0064] Fig. 6C is a simplified flowchart illustrating operation of a 
dicer-cut location detector constructed and operative in 
accordance with a preferred embodiment of the present 
invention; 

[0065] Fig. 7 A is a simplified block diagram of a target gene 

binding-site detector constructed and operative in accor- 
dance with a preferred embodiment of the present inven- 
tion; 

[0066] Fig. 7B is a simplified flowchart illustrating operation of a 



target gene binding-site detector constructed and opera- 
tive in accordance with a preferred embodiment of the 
present invention; 
[0067] pig. 8 is a simplified flowchart illustrating operation of a 
function & utility analyzer constructed and operative in 
accordance with a preferred embodiment of the present 
invention; 

[0068] pig. 9 is a simplified diagram describing a novel bioinfor- 
matically detected group of regulatory polynucleotides, 
referred to here as Genomic Record (GR) polynucleotides, 
each of which encodes an x operon-like x cluster of novel 
miRNA-like oligonucleotides, which in turn modulate ex- 
pression of one or more target genes; 

[0069] Fig. 10 is a block diagram illustrating different utilities of 
novel oligonucleotides and novel operon-like polynu- 
cleotides, both of the present invention; 

[0070] Figs. 11A and 11B are simplified diagrams, which when 
taken together illustrate a mode of oligonucleotide ther- 
apy applicable to novel oligonucleotides of the present in- 
vention; 

[0071] Fig. 12A is a histogram representing the distribution of 
known miRNA oligonucleotides and that of predicted 
miRNA-like hairpin structures extracted from expressed 



genome sequences with respect to their hairpin detector 
score; 

[0072] pig. 12B is a table summarizing laboratory validation re- 
sults which validate efficacy of a bioinformatic oligonu- 
cleotide detection system constructed and operative in 
accordance with a preferred embodiment of the present 
invention; 

[0073] pig. 13A and Fig. 13B are a picture and a summary table 
of laboratory results validating the expression of novel 
human oligonucleotides detected by a bioinformatic 
oligonucleotide detection engine constructed and opera- 
tive in accordance with a preferred embodiment of the 
present invention, thereby validating the efficacy of the 
oligonucleotide detection engine of the present invention; 

[0074] Fig. 14A is a schematic representation of an "operon-like" 
cluster of novel human hairpin sequences detected bioin- 
formatically by a bioinformatic oligonucleotide detection 
engine constructed and operative in accordance with a 
preferred embodiment of the present invention, and non- 
GAM hairpin useful as negative controls thereto; 

[0075] Fig. 14B is a schematic representation of secondary fold- 
ing of hairpins of the operon-like cluster of Fig. 14A; 

[0076] Fig. 14C is a picture of laboratory results demonstrating 



expression of novel oligonucleotides of Figs. 14A and 
14B, and lack of expression of the negative controls, 
thereby validating efficacy of bioinformatic detection of 
GAM oligonucleotides and GR polynucleotides detected by 
a bioinformatic oligonucleotide detection engine, con- 
structed and operative in accordance with a preferred em- 
bodiment of the present invention; 
[0077] pig. 15A is an annotated sequence of EST72223 compris- 
ing known human miRNA oligonucleotide MIR98 and novel 
human oligonucleotide GAM25 PRECURSOR detected by 
the oligonucleotide detection system of the present in- 
vention; and 

[0078] pigs. 15B, 15C and 15D are pictures of laboratory results 
demonstrating laboratory confirmation of expression of 
known human oligonucleotide MIR98 and of novel bioin- 
formatically detected human GAM25 RNA respectively, 
both of Fig. 15A, thus validating the bioinformatic 
oligonucleotide detection system of the present invention; 

[0079] Fig. 16 presents pictures of laboratory results demon- 
strating laboratory confirmation of 'dicing' of four novel 
bioinformatically detected HIV1 GAM PRECURSORSs into 
their corresponding mature GAM RNAs, herein designated 
VGAM2032-B (Fig. 16B),VGAM3249-A (Fig. 16C), 



GAM507-B (Fig. 16D) and VGAM1016-B (Fig.l6E); 
Brief Description of Sequences 



[0080] a Sequence Listing of genomic sequences of the present 
invention designated SEQ ID NO:l through SEQ ID: 
399737 is attached to this application, and is hereby in- 
corporated herein. The genomic listing comprises the fol- 
lowing nucleotide sequences: nucleotide sequences of 
15484 GAM precursors of respective novel oligonu- 
cleotides of the present invention; nucleotide sequences 
of 14452 GAM RNA oligonucleotides of respective novel 
oligonucleotides of the present invention; and nucleotide 
sequences of 127998 target gene binding sites of respec- 
tive novel oligonucleotides of the present invention. 
Detailed Description 

[0081] Reference is now made to Fig. 1, which is a simplified dia- 
gram describing a plurality of novel bioinformatically de- 
tected oligonucleotide of the present invention referred to 
here as Genomic Address Messenger (GAM) oligonu- 
cleotide, which modulates expression of respective target 
genes whose function and utility is known in the art. 

[0082] gam is a novel bioinformatically detectable regulatory, 
non protein coding, microRNA (miRNA)-like oligonu- 



cleotide. The method by which CAM is detected is de- 
scribed with additional reference to Figs. 1-8. 

[0083] GAM PRECURSOR is preferably encoded by a viral genome. 
Alternatively or additionally GAM PRECURSOR is encoded 
by the human genome. GAM TARGET GENE is a gene en- 
coded by the human genome. Alternatively or additionally 
GAM TARGET GENE is a gene encoded by a viral genome. 

[0084] GAM PRECURSOR encodes a GAM PRECURSOR RNA. Similar 
to other miRNA oligonucleotides GAM PRECURSOR RNA 
does not encode a protein. 

[0085] GAM PRECURSOR RNA folds onto itself, forming GAM 
FOLDED PRECURSOR RNA which has a two-dimensional 
N hairpin structure \ As is well known in the art, this 
" hairpin structure \ is typical of RNA encoded by known 
miRNA oligonucleotides, and is due to the fact that the 
nucleotide sequence of the first half of the RNA encoded 
by a miRNA oligonucleotide is a fully or partially comple- 
mentary sequence of the nucleotide sequence of the sec- 
ond half thereof. By complementary is meant a sequence 
which is reversed and wherein each nucleotide is replaced 
by a complementary nucleotide, as is well known in the art 
(e.g. ATGGC is the complementary sequence of GCCAT). 
An enzyme complex designated DICER COMPLEX, x dices x 



the GAM FOLDED PRECURSOR RNA yielding a GAM RNA in 
the form of a single stranded ~22nt long RNA segment. As 
is known in the art, ^dicing^ of a hairpin structured RNA 
precursor product to yield a short ~22nt RNA segment is 
catalyzed by an enzyme complex Dicer RNaselll together 
with other necessary proteins. 
[0086] gam TARGET GENE encodes a corresponding messenger 

RNA, GAM TARGET RNA. Each of the GAM TARGET RNAs of 
the present invention comprises three regions, as is typi- 
cal of mRNA of a protein coding gene: a 5 x untranslated 
region, a protein coding region and a 3 X untranslated re- 
gion, designated 5 V UTR, PROTEIN CODING and 3^UTR re- 
spectively. 

[0087] gam RNA binds complementarily to one or more target 
binding sites located in untranslated regions of each of 
the GAM TARGET RNAs of the present invention. This 
complementary binding is due to the fact that the nu- 
cleotide sequence of GAM RNA is a partial or fully comple- 
mentary sequence of the nucleotide sequence of each of 
the target binding sites. As an illustration, Fig. 1 shows 
three such target binding sites, designated BINDING SITE I, 
BINDING SITE II and BINDING SITE III respectively. It is ap- 
preciated that the number of target binding sites shown in 



Fig.l is only illustrative and that any suitable number of 
target binding sites may be present. It is further appreci- 
ated that although Fig.l shows target binding sites only in 
the 3 X UTR region these target binding sites may be lo- 
cated instead in the 5 x UTR region or in both 3 x UTR and 
5 X UTR regions. 

[0088] The complementary binding of GAM RNA to target binding 
sites on GAM TARGET RNA, such as BINDING SITE I, BIND- 
ING SITE II and BINDING SITE III, inhibits translation of 
each of the GAM TARGET RNAs of the present invention 
into GAM TARGET PROTEIN, which is shown surrounded 
by a broken line. 

[0089] | t j S appreciated that GAM TARGET GENE in fact represents 
a plurality of GAM target genes. The mRNA of each one of 
this plurality of GAM target genes comprises one or more 
target binding sites, each having a nucleotide sequence 
which is at least partly complementary to GAM RNA and 
which when bound by GAM RNA causes inhibition of 
translation of the GAM target mRNA into a corresponding 
GAM target protein. 

[0090] The mechanism of the translational inhibition exerted by 
GAM RNA on one or more GAM TARGET GENEs may be 
similar or identical to the known mechanism of transla- 



tional inhibition exerted by known miRNA oligonu- 
cleotides. 

[0091] Nucleotide sequences of each of a plurality of GAM 

oligonucleotides described by Fig. 1 and their respective 
genomic sources and genomic locations are set forth in 
Tables 1-3, hereby incorporated herein. 

[0092] Nucleotide sequences of CAM PRECURSOR RNAs, and a 

schematic representation of a predicted secondary folding 
of GAM FOLDED PRECURSOR RNAs, of each of a plurality 
of GAM oligonucleotides described by Fig. 1 are set forth 
in Table 4 , hereby incorporated herein. 

[0093] Nucleotide sequences of diced x GAM RNAs of each of a 
plurality of GAM oligonucleotides described by Fig. 1 are 
set forth in Table 5, hereby incorporated herein. 

[0094] Nucleotide sequences of target binding sites, such as 

BINDING SITE I, BINDING SITE II and BINDING SITE III found 
on GAM TARGET RNAs of each of a plurality of GAM 
oligonucleotides described by Fig. 1, and a schematic rep- 
resentation of the complementarity of each of these target 
binding sites to each of a plurality of GAM RNAs described 
by Fig. 1 are set forth in Tables 6-7, hereby incorporated 
herein. 

[0095] it is appreciated that specific functions and accordingly 



utilities of each of a plurality of GAM oligonucleotides de- 
scribed by Fig. 1 are correlated with, and may be deduced 
from the identity of the GAM TARGET GENEs inhibited 
thereby, and whose functions are set forth in Table 8, 
hereby incorporated herein. 

[0096] studies documenting well known correlations between 
each of a plurality of GAM TARGET GENEs of the GAM 
oligonucleotides of Fig. 1, and known functions and dis- 
eases are listed in Table 9, hereby incorporated herein. 

[0097] The present invention discloses a novel group of viral and 
human oligonucleotides, belonging to the miRNA-like 
oligonucleotides group, here termed GAM oligonu- 
cleotides, for which a specific complementary binding has 
been determined bioinformatically. 

[0098] Reference is now made to Fig. 2 which is a simplified 

block diagram illustrating a bioinformatic oligonucleotide 
detection system and method constructed and operative 
in accordance with a preferred embodiment of the present 
invention. 

[0099] An important feature of the present invention is a BIOIN- 
FORMATIC OLIGONUCLEOTIDE DETECTION ENGINE 100, 
which is capable of bioinformatically detecting oligonu- 
cleotides of the present invention. 



[0100] The functionality of the BIOINFORMATIC OLIGONU- 
CLEOTIDE DETECTION ENGINE 100 includes receiving EX- 
PRESSED RNA DATA 102, SEQUENCED DNA DATA 104, and 
PROTEIN FUNCTION DATA 106; performing a complex 
process of analysis of this data as elaborated hereinbelow, 
and based on this analysis provides information, desig- 
nated by reference numeral 108, identifying and describ- 
ing features of novel oligonucleotides. 

[0101] EXPRESSED RNA DATA 102 comprises published ex- 
pressed sequence tags (EST) data, published mRNA data, 
as well as other published RNA data. SEQUENCED DNA 
DATA 104 comprises alphanumeric data representing ge- 
nomic sequences and preferably including annotations 
such as information indicating the location of known pro- 
tein coding regions relative to the genomic sequences. 

[0102] PROTEIN FUNCTION DATA 106 comprises information 

from scientific publications e.g. physiological functions of 
known proteins and their connection, involvement and 
possible utility in treatment and diagnosis of various dis- 
eases. 

[0103] EXPRESSED RNA DATA 102 and SEQUENCED DNA DATA 
104 may preferably be obtained from data published by 
the National Center for Biotechnology Informatiion (NCBI) 



at the National Institute of Health (NIH) Qenuth J. P. (2000). 
Methods Mol. Biol. 132:301-312(2000), herein incorpo- 
rated by reference) as well as from various other pub- 
lished data sources. PROTEIN FUNCTION DATA 106 may 
preferably be obtained from any one of numerous relevant 
published data sources, such as the Online Mendelian In- 
herited Disease In Man (OMIM(TM), Hamosh et al., Nucleic 
Acids Res. 30: 52-55(2002)) database developed by John 
Hopkins University, and also published by NCBI (2000). 

[0104] p r j 0r t0 or during actual detection of BIOINFORMATICALLY 
DETECTED GROUP OF NOVEL OLIGONUCLEOTIDES 108 by 
the BIOINFORMATIC OLIGONUCLEOTIDE DETECTION EN- 
GINE 100, BIOINFORMATIC OLIGONUCLEOTIDE DETECTION 
ENGINE TRAINING & VALIDATION FUNCTIONALITY 110 is 
operative. This functionality uses one or more known 
miRNA oligonucleotides as a training set to train the 
BIOINFORMATIC OLIGONUCLEOTIDE DETECTION ENGINE 
100 to bioinformatically recognize miRNA-like oligonu- 
cleotides, and their respective potential target binding 
sites. BIOINFORMATIC OLIGONUCLEOTIDE DETECTION EN- 
GINE TRAINING & VALIDATION FUNCTIONALITY 110 is fur- 
ther described hereinbelow with reference to Fig. 3. 

[° 1 05] The BIOINFORMATIC OLIGONUCLEOTIDE DETECTION EN- 



GINE 100 preferably comprises several modules which are 
preferably activated sequentially, and are described as 
follows: 

[0 1 06] a NON-CODING GENOMIC SEQUENCE DETECTOR 112 op- 
erative to bioinformatically detect non-protein coding ge- 
nomic sequences. The NON-CODING GENOMIC SEQUENCE 
DETECTOR 112 is further described herein below with ref- 
erence to Figs. 4A and 4B. 

[0 1 ° 7 ] A HAIRPIN DETECTOR 114 operative to bioinformatically 
detect genomic x hairpin-shaped x sequences, similar to 
GAM FOLDED PRECURSOR RNA (Fig. 1). The HAIRPIN DE- 
TECTOR 114 is further described herein below with refer- 
ence to Figs. 5A and 5B. 

[0 1 °8] a DICER-CUT LOCATION DETECTOR 116 operative to 
bioinformatically detect the location on a GAM FOLDED 
PRECURSOR RNA which is enzymatically cut by DICER 
COMPLEX (Fig. 1), yielding diced GAM RNA. The DICER- 
CUT LOCATION DETECTOR 116 is further described herein 
below with reference to Figs. 6A - 6C. 

[° 1 09] A TARGET GENE BINDING-SITE DETECTOR 118 operative 
to bioinformatically detect target genes having binding 
sites, the nucleotide sequence of which is partially com- 
plementary to that of a given genomic sequence, such as 



a nucleotide sequence cut by DICER COMPLEX. The TAR- 
GET GENE BINDING-SITE DETECTOR 118 is further de- 
scribed hereinbelow with reference to Figs. 7 A and 7B. 
[0 1 1 0] A FUNCTION & UTILITY ANALYZER 120 operative to ana- 
lyze the function and utility of target genes in order to 
identify target genes which have a significant clinical 
function and utility. The FUNCTION & UTILITY ANALYZER 
120 is further described hereinbelow with reference to 
Fig. 8. 

[° 1 1 1 ] According to a preferred embodiment of the present in- 
vention the engine 100 may employ a cluster of 40 PCs ( 
XEON (R) , 2.8GHz, with 80GB storage each), connected by 
Ethernet to 8 servers (2-CPU, XEON (TM) 1.2-2.2GHz, with 
~200GB storage each), combined with an 8-processor 
server (8-CPU, Xeon 550Mhz w/ 8GB RAM) connected via 
2 HBA fiber-channels to an EMC CLARIION (TM) 
100-disks, 3.6 Terabyte storage device. A preferred em- 
bodiment of the present invention may also preferably 
comprise software which utilizes a commercial database 
software program, such as MICROSOFT (TM) SQL Server 
2000. It is appreciated that the above mentioned hard- 
ware configuration is not meant to be limiting, and is 
given as an illustration only. The present invention may be 



implemented in a wide variety of hardware and software 
configurations. 

[0112] The present invention discloses 15484 novel oligonu- 
cleotides of the GAM group of oligonucleotides, which 
have been detected bioinformatically and 699 novel 
polynucleotides of the GR group of polynucleotides , 
which have been detected bioinformatically. Laboratory 
confirmation of bioinformatically predicted oligonu- 
cleotides of the GAM group of oligonucleotides, and sev- 
eral bioinformatically predicted polynucleotides of the GR 
group of polynucleotides, is described hereinbelow with 
reference to Figs. 12 - 15D. Laboratory confirmation of 
bioinformatically predicted oligonucleotides of the viral 
GAM group of oligonucleotides , and several bioinformati- 
cally predicted viral polynucleotides of the GR group of 
polynucleotides, is described hereinbelow with reference 
to Fig. 16. 

[0113] Reference is now made to Fig. 3 which is a simplified 

flowchart illustrating operation of a preferred embodiment 
of the BIOINFORMATIC OLIGONUCLEOTIDE DETECTION 
ENGINE TRAINING &VALIDATION FUNCTIONALITY 110 de- 
scribed hereinabove with reference to Fig. 2. 

[O 114 ] BIOINFORMATIC OLIGONUCLEOTIDE DETECTION ENGINE 



TRAINING & VALIDATION FUNCTIONALITY 110 begins by 
training the BIOINFORMATIC OLIGONUCLEOTIDE DETEC- 
TION ENGINE 100 (Fig. 2) to recognize one or more known 
miRNA oligonucleotides, as designated by reference nu- 
meral 122. This training step comprises HAIRPIN DETEC- 
TOR TRAINING & VALIDATION FUNCTIONALITY 124, fur- 
ther described hereinbelow with reference to Fig. 5A, 
DICER-CUT LOCATION DETECTOR TRAINING & VALIDA- 
TION FUNCTIONALITY 126, further described hereinbelow 
with reference to Fig. 6A and 6B, and TARGET GENE BIND- 
ING-SITE DETECTOR TRAINING & VALIDATION FUNCTION- 
ALITY 128, further described hereinbelow with reference 
to Fig. 7A. 

15 ] Next, the BIOINFORMATIC OLIGONUCLEOTIDE DETECTION 
ENGINE TRAINING & VALIDATION FUNCITONALITY 110 is 
operative bioinformatically detect novel oligonucleotides, 
using BIOINFORMATIC OLIGONUCLEOTIDE DETECTION EN- 
GINE 100 (Fig. 2), as designated by reference numeral 
130. Wet lab experiments are preferably conducted in or- 
der to validate expression and preferably function of some 
samples of the novel oligonucleotides detected by the 
BIOINFORMATIC OLIGONUCLEOTIDE DETECTION ENGINE 
100, as designated by reference numeral 132. Figs. 13A - 



15D illustrate examples of wet-lab validation of sample 
novel human oligonucleotides bioinformatically detected 
in accordance with a preferred embodiment of the present 
invention. Laboratory confirmation of bioinformatically 
predicted oligonucleotides of the viral GAM group of 
oligonucleotides , and several bioinformatically predicted 
viral polynucleotides of the GR group of polynucleotides, 
is described hereinbelow with reference to Fig. 16. 
16 ] Reference is now made to Fig. 4A which is a simplified 
block diagram of a preferred implementation of the NON- 
CODING GENOMIC SEQUENCE DETECTOR 112 described 
hereinabove with reference to Fig. 2. The NON-PROTEIN 
CODING GENOMIC SEQUENCE DETECTOR 112 preferably 
receives at least two types of published genomic data: EX- 
PRESSED RNA DATA 102 and SEQUENCED DNA DATA 104. 
The EXPRESSED RNA DATA 102 may include, inter alia, EST 
data, EST clusters data, EST genome alignment data and 
mRNA data. Sources for EXPRESSED RNA DATA 102 in- 
clude NCBI dbEST, NCBI UniGene clusters and mapping 
data, and TIGR gene indices (Kirkness F. and Kerlavage, 
A.R., Methods Mol. Biol. 69:261-268 (1997)). SEQUENCED 
DNA DATA 104 may include sequence data (FASTA format 
files), and feature annotations (GenBank file format) 



mainly from NCBI databases. Based on the above men- 
tioned input data, the NON-PROTEIN CODING GENOMIC 
SEQUENCE DETECTOR 112 produces a plurality of NON- 
PROTEIN CODING GENOMIC SEQUENCES 136. Preferred 
operation of the NON-PROTEIN CODING GENOMIC SE- 
QUENCE DETECTOR 112 is described hereinbelow with 
reference to Fig. 4B. 
[0117] Reference is now made to Fig. 4B which is a simplified 
flowchart illustrating a preferred operation of the NON- 
CODING GENOMIC SEQUENCE DETECTOR 112 of Fig. 2. 
Detection of NON-PROTEIN CODING GENOMIC SEQUENCES 
136, generally preferably progresses along one of the fol- 
lowing two paths: 

A first path for detecting NON-PROTEIN CODING GENOMIC 
SEQUENCES 136 (Fig. 4A) begins with receipt of a plurality 
of known RNA sequences, such as EST data. Each RNA se- 
quence is first compared with known protein-coding DNA 
sequences, in order to select only those RNA sequences 
which are non-protein coding, i.e. intergenic or intronic 
sequences. This can preferably be performed by using one 
of many alignment algorithms known in the art, such as 
BLAST (Altschul et al.J. Mol. Biol. 215:403-410 (1990)). 
This sequence comparison preferably also provides local- 



ization of the RNA sequence on the DNA sequences. 

[0119] Alternatively, selection of non-protein coding RNA se- 
quences and their localization on the DNA sequences can 
be performed by using publicly available EST cluster data 
and genomic mapping databases, such as the UNIGENE 
database published by NCBI or the TIGR database. Such 
databases, map expressed RNA sequences to DNA se- 
quences encoding them, find the correct orientation of 
EST sequences, and indicate mapping of ESTs to protein 
coding DNA regions, as is well known in the art. Public 
databases, such as TIGR, may also be used to map an EST 
to a cluster of ESTs, known in the art as Tentative Human 
Consensus and assumed to be expressed as one segment. 
Publicly available genome annotation databases, such as 
NCBIs GenBank, may also be used to deduce expressed 
intronic sequences. 

[0120] Optionally, an attempt may be made to "expand" the non- 
protein RNA sequences thus found, by searching for tran- 
scription start and end signals, respectively upstream and 
downstream of the location of the RNA on the DNA, as is 
well known in the art. 

[° 121 ] A second path for detecting NON-PROTEIN CODING GE- 
NOMIC SEQUENCES 136 (Fig. 4A) begins with receipt of 



DNA sequences. The DNA sequences are parsed into non 
protein coding sequences, using published DNA annota- 
tion data, by extracting those DNA sequences which are 
between known protein coding sequences. Next, tran- 
scription start and end signals are sought. If such signals 
are found, and depending on their robustness, probable 
expressed non-protein coding genomic sequences are 
obtained. Such approach is especially useful for identify- 
ing novel GAM oligonucleotides which are found in prox- 
imity to other known miRNA oligonucleotides, or other 
wet-lab validated GAM oligonucleotides. Since, as de- 
scribed hereinbelow with reference to Fig. 9, GAM 
oligonucleotides are frequently found in clusters; se- 
quences located near known miRNA oligonucleotides are 
more likely to contain novel GAM oligonucleotides. Op- 
tionally, sequence orthology, i.e. sequence conservation in 
an evolutionary related species, may be used to select ge- 
nomic sequences having a relatively high probability of 
containing expressed novel GAM oligonucleotides. It is 
appreciated that in detecting non-human GAM oligonu- 
cleotides of the present invention the BIOINFORMATIC 
OLIGONUCLEOTIDE DETECTION ENGINE 100 utilizes the 
input genomic sequences, without filtering protein coding 



regions detected by the NON-CODING GENOMIC SE- 
QUENCE DETECTOR 112, hence NON-PROTEIN CODING 
GENOMIC SEQUENCES 136 referes to GENOMIC SE- 
QUENCES only. 

[0122] Reference is now made to Fig. 5A which is a simplified 

block diagram of a preferred implementation of the HAIR- 
PIN DETECTOR 114 described hereinabove with reference 
to Fig. 2. 

[0 12 3] The goal of the HAIRPIN DETECTOR 114 is to detect hair- 
pin-shaped genomic sequences, similar to those of known 
miRNA oligonucleotides. A hairpin- shaped genomic se- 
quence is a genomic sequence, having a first half which is 
at least partially complementary to a second half thereof, 
which causes the halves to folds onto themselves, thereby 
forming a hairpin structure, as mentioned hereinabove 
with reference to Fig. 1. 

[O 124 ] The HAIRPIN DETECTOR 114 (Fig. 2) receives a plurality of 
NON-PROTEIN CODING GENOMIC SEQUENCES 136 (Fig. 
4A). Following operation of HAIRPIN DETECTOR TRAINING 
& VALIDATION FUNCTIONALITY 124 ( Fig. 3), the HAIRPIN 
DETECTOR 114 is operative to detect and output hairpin- 
shaped sequences, which are found in the NON-PROTEIN 
CODING GENOMIC SEQUENCES 136. The hairpin-shaped 



sequences detected by the HAIRPIN DETECTOR 114 are 
designated HAIRPINS STRUCTURES ON GENOMIC SE- 
QUENCES 138. A preferred mode of operation of the 
HAIRPIN DETECTOR 114 is described hereinbelow with 
reference to Fig. 5B. 
[0 12 5] HAIRPIN DETECTOR TRAINING & VALIDATION FUNCTION- 
ALITY 124 includes an iterative process of applying the 
HAIRPIN DETECTOR 114 to known hairpin shaped miRNA 
precursor sequences, calibrating the HAIRPIN DETECTOR 
114 such that it identifies a training set of known hairpin- 
shaped miRNA precursor sequences, as well as other sim- 
ilarly hairpin-shaped sequences. In a preferred embodi- 
ment of the present invention, the HAIRPIN DETECTOR 
TRAINING & VALIDATION FUNCTIONALITY 124 trains the 
HAIRPIN DETECTOR 114 and validates each of the steps of 
operation thereof described hereinbelow with reference to 
Fig. 5B 

[0 12 6] The HAIRPIN DETECTOR TRAINING & VALIDATION FUNC- 
TIONALITY 124 preferably uses two sets of data: the 
aforesaid training set of known hairpin-shaped miRNA 
precursor sequences, such as hairpin-shaped miRNA pre- 
cursor sequences of 440 miRNA oligonucleotides of H. 
sapiens, M. musculus, C. elegans, C. Brigssae and D. 



Melanogaster, annotated in the RFAM database 
(Griffiths-Jones 2003), and a background set of about 
1000 hairpin-shaped sequences found in expressed non- 
protein coding human genomic sequences. The back- 
ground set is expected to comprise some valid, previously 
undetected hairpin-shaped miRNA-like precursor se- 
quences, and many hairpin-shaped sequences which are 
not hairpin-shaped miRNA-like precursors. 

[° 127 ] In a preferred embodiment of the present invention the 
efficacy of the HAIRPIN DETECTOR 114 (Fig. 2) is con- 
firmed. For example, when a similarity threshold is chosen 
such that 87% of the known hairpin-shaped miRNA pre- 
cursors are successfully predicted, only 21.8% of the 1000 
background set of hairpin-shaped sequences are pre- 
dicted to be hairpin-shaped miRNA-like precursors. 

[0128] Reference is now made to Fig. 5B which is a simplified 
flowchart illustrating preferred operation of the HAIRPIN 
DETECTOR 114 of Fig. 2. The HAIRPIN DETECTOR 114 
preferably initially uses a secondary structure folding al- 
gorithm based on free-energy minimization, such as the 
MFOLD algorithm, described in Mathews et al. J. Mol. Biol. 
288:911-940 (1999) and Zuker, M. Nucleic Acids Res. 31: 
3406-3415 (2003), the disclosure of which is hereby in- 



corporated by reference. This algorithm is operative to 
calculate probable secondary structure folding patterns of 
the NON-PROTEIN CODING GENOMIC SEQUENCES 136 
(Fig. 4A) as well as the free-energy of each of these prob- 
able secondary folding patterns. The secondary structure 
folding algorithm, such as the MFOLD algorithm 
(Mathews, 1997; Zuker 2003), typically provides a listing 
of the base-pairing of the folded shape, i.e. a listing of 
each pair of connected nucleotides in the sequence. 
[0 12 9] Next, the HAIRPIN DETECTOR 114 analyzes the results of 
the secondary structure folding patterns, in order to de- 
termine the presence and location of hairpin folding 
structures. The goal of this second step is to assess the 
base-pairing listing provided by the secondary structure 
folding algorithm, in order to determine whether the 
base-pairing listing describes one or more hairpin type 
bonding pattern. Preferably, sequence segment corre- 
sponding to a hairpin structure is then separately ana- 
lyzed by the secondary structure folding algorithm in or- 
der to determine its exact folding pattern and free-en- 
ergy. 

[0 1 30] The HAIRPIN DETECTOR 114 then assesses the hairpin 

structures found by the previous step, comparing them to 



hairpin structures of known miRNA precursors, using vari- 
ous characteristic hairpin structure features such as its 
free-energy and its thermodynamic stability, the amount 
and type of mismatched nucleotides and the existence of 
sequence repeat-elements, number of mismatched nu- 
cleotides in positions 18-22 counting from loop, and Per- 
cent of G nucleotide. Only hairpins that bear statistically 
significant resemblance to the training set of hairpin 
structures of known miRNA precursors, according to the 
abovementioned parameters, are accepted. 
[° 131 ] In a preferred embodiment of the present invention, simi- 
larity to the training set of hairpin structures of known 
miRNA precursors is determined using a "similarity score" 
which is calculated using a multiplicity of terms, where 
each term is a function of one of the abovementioned 
hairpin structure features. The parameters of each func- 
tion are found heuristically from the set of hairpin struc- 
tures of known miRNA precursors, as described herein- 
above with reference to HAIRPIN DETECTOR TRAINING & 
VALIDATION FUNCTIONALITY 124 (Fig. 3). The selection of 
the features and their function parameters is optimized so 
as to achieve maximized separation between the distribu- 
tion of similarity scores validated miRNA-precursor hair- 



pin structures, and the distribution of similarity scores of 
hairpin structures detected in the background set men- 
tioned hereinabove with reference to Fig. 5B. 

[° 132 ] In an alternative preferred embodiment of the present in- 
vention, the step described in the preceding paragraph 
may be split into two stages. A first stage implements a 
simplified scoring method, typically based on threshold- 
ing a subset of the hairpin structure features described 
hereinabove, and may employ a minimum threshold for 
hairpin structure length and a maximum threshold for 
free energy. A second stage is preferably more stringent, 
and preferably employs a full calculation of the weighted 
sum of terms described hereinabove. The second stage 
preferably is performed only on the subset of hairpin 
structures that survived the first stage. 

[° 1 33] The HAIRPIN DETECTOR 114 also attempts to select hair- 
pin structures whose thermodynamic stability is similar to 
that of hairpin structures of known miRNA precursors. 
This may be achieved in various ways. A preferred em- 
bodiment of the present invention utilizes the following 
methodology, preferably comprising three logical steps: 

[° 134 ] First, the HAIRPIN DETECTOR 114 attempts to group hair- 
pin structures into "families" of closely related hairpin 



structures. As is known in the art, a secondary structure 
folding algorithm typically provides multiple alternative 
folding patterns, for a given genomic sequence and indi- 
cates the free energy of each alternative folding pattern.. 
It is a particular feature of the present invention that the 
HAIRPIN DETECTOR 114 preferably assesses the various 
hairpin structures appearing in the various alternative 
folding patterns and groups x hairpin structures which ap- 
pear at identical or similar sequence locations in various 
alternative folding patterns into common sequence loca- 
tion based "families" of hairpins. For example, all hairpin 
structures whose center is within 7 nucleotides of each 
other may be grouped into a family". Hairpin structures 
may also be grouped into a family" if their nucleotide se- 
quences are identical or overlap to a predetermined de- 
gree. 

[0135] ^ is also a particular feature of the present invention that 
the hairpin structure "families" are assessed in order to 
select only those families which represent hairpin struc- 
tures that are as thermodynamically stable as those of 
hairpin structures of known miRNA precursors. Preferably 
only families which are represented in at least a selected 
majority of the alternative secondary structure folding 



patterns, typically 65%, 80% or 100% are considered to be 
sufficiently stable. Our tests suggest that only about 50% 
of the hairpin structures, predicted by the MFOLD algo- 
rithm with default parameters, are members of sufficiently 
stable families, comparing to about 90% of the hairpin 
structures that contain known miRNAs. This percent de- 
pends on the size of the fraction that was fold. In an alter- 
native embodiment of the present invention we use frac- 
tions of size lOOOnts as preferable size. Different embod- 
iment uses other sizes of genomics sequences, more or 
less strict demand for representation in the alternative 
secondary structure folding patterns. 
[0136] it is an additional particular feature of the present inven- 
tion that the most suitable hairpin structure is selected 
from each selected family. For example, a hairpin struc- 
ture which has the greatest similarity to the hairpin struc- 
tures appearing in alternative folding patterns of the fam- 
ily may be preferred. Alternatively or additionally, the 
hairpin structures having relatively low free energy may be 
preferred. 

[0137] Alternatively or additionally considerations of homology to 
hairpin structures of other organisms and the existence of 
clusters of thermodynamically stable hairpin structures 



located adjacent to each other along a sequence may be 
important in selection of hairpin structures. The tightness 
of the clusters in terms of their location and the occur- 
rence of both homology and clusters may be of signifi- 
cance. 

[0138] Reference is now made to Figs. 6A - 6C which together 
describe the structure and operation of the DICER-CUT 
LOCATION DETECTOR 116, described hereinabove with 
Fig. 2. 

[0139] pig. 6A is a simplified block diagram of a preferred imple- 
mentation of the DICER-CUT LOCATION DETECTOR 116. 
The goal of the DICER-CUT LOCATION DETECTOR 116 is 
to detect the location in which the DICER COMPLEX, de- 
scribed hereinabove with reference to Fig. 1, dices CAM 
FOLDED PRECURSOR RNA, yielding CAM RNA. 

[0 14 0] The DICER-CUT LOCATION DETECTOR 116 therefore re- 
ceives a plurality of HAIRPIN STRUCTURES ON GENOMIC 
SEQUENCES 138 (Fig. 5A), and following operation of 
DICER-CUT LOCATION DETECTOR TRAINING & VALIDA- 
TION FUNCTIONALITY 126 (Fig. 3), is operative to detect a 
plurality of DICER-CUT SEQUENCES FROM HAIRPIN STRUC- 
TURES 140. 

[0141] Reference is now made to Fig. 6B which is a simplified 



flowchart illustrating a preferred implementation of 
DICER-CUT LOCATION DETECTOR TRAINING & VALIDA- 
TION FUNCTIONALITY 126. 
[0 142 ] A general goal of the DICER-CUT LOCATION DETECTOR 
TRAINING & VALIDATION FUNCTIONALITY 126 is to ana- 
lyze the dicer-cut locations of known diced miRNA on re- 
spective hairpin shaped miRNA precursors in order to de- 
termine a common pattern in these locations, which can 
be used to predict dicer cut locations on GAM folded pre- 
cursor RNAs. 

[0143] The dicer-cut locations of known miRNA precursors are 
obtained and studied. Locations of the 5' and/or 3' ends 
of the known diced miRNAs are preferably represented by 
their respective distances from the 5' end of the corre- 
sponding hairpin shaped miRNA precursor. Additionally or 
alternatively, the 5' and/or 3' ends of the known diced 
miRNAs are preferably represented by the relationship be- 
tween their locations and the locations of one or more nu- 
cleotides along the hairpin shaped miRNA precursor. Ad- 
ditionally or alternatively, the 5' and/or 3' ends of the 
known diced miRNAs are preferably represented by the 
relationship between their locations and the locations of 
one or more bound nucleotide pairs along the hairpin 



shaped miRNA precursor. Additionally or alternatively, the 
5' and/or 3' ends of the known diced miRNAs are prefer- 
ably represented by the relationship between their loca- 
tions and the locations of one or more mismatched nu- 
cleotide pairs along the hairpin shaped miRNA precursor. 
Additionally or alternatively, the 5' and/or 3' ends of the 
known diced miRNAs are preferably represented by the 
relationship between their locations and the locations of 
one or more unmatched nucleotides along the hairpin 
shaped miRNA precursor. Additionally or alternatively, lo- 
cations of the 5' and/or 3' ends of the known diced miR- 
NAs are preferably represented by their respective dis- 
tances from the loop located at the center of the corre- 
sponding hairpin shaped miRNA precursor. 

[0144] one or more of the foregoing location metrics may be 

employed in the training and validation functionality. Ad- 
ditionally, metrics related to the nucleotide content of the 
diced miRNA and/or of the hairpin shaped miRNA precur- 
sor may be employed. 

[° 145 ] In a preferred embodiment of the present invention, 

DICER-CUT LOCATION DETECTOR TRAINING & VALIDA- 
TION FUNCTIONALITY 126 preferably employs standard 
machine learning techniques known in the art of machine 



learning for analysis of existing patterns in a given "train- 
ing set" of examples. These techniques are capable, to a 
certain degree, of detecting similar patterns in other, pre- 
viously unseen examples. Such machine learning tech- 
niques include, but are not limited to neural networks, 
Bayesian networks, Support Vector Machines (SVM), Ge- 
netic Algorithms, Markovian modeling, Maximum Likeli- 
hood modeling, Nearest Neighbor algorithms, Decision 
trees and other techniques, as is well known in the art. 
[0146] | n accordance with one embodiment of the present inven- 
tion, machine learning predictors, such as a Support Vec- 
tor Machine (SVM) predictor, are applied to the aforemen- 
tioned training set and are operative, for example to test 
every possible nucleotide on a hairpin as a candidate for 
being the 5' end or the 3' end of a diced CAM RNA. More 
preferred machine learning predictors include predictors 
based on Nearest Neighbor, Bayesian modeling, and K- 
nearest-neighbor algorithms. A training set of the known 
miRNA precursor sequences is preferably used for training 
multiple separate classifiers or predictors, each of which 
produces a model for the 5' and/or 3' end locations of a 
diced miRNA with respect to its hairpin precursor. The 
models take into account one or more of the various 



miRNA location metrics described above. 

[0147] Performance of the resulting predictors, evaluated on the 
abovementioned validation set of 440 published miRNAs 
using k-fold cross validation (Mitchell, 1997) with k = 3, 
is found to be as follows: in 70% of known miRNAs 5'-end 
location is correctly determined by an SVM predictor 
within up to 2 nucleotides; a Nearest Neighbor (EDIT DIS- 
TANCE) predictor achieves 56% accuracy (247/440); a 
Two-Phased predictor that uses Bayesian modeling (TWO 
PHASED) achieves 80% accuracy (352/440), when only the 
first phase is used. When the second phase (strand choice) 
is implemented by a nave Bayesian model the accuracy is 
55% (244/440), and when the K-nearest-neighbor model- 
ing is used for the second phase, 374/440 decision are 
made and the accuracy is 65% (242/374). A K-near- 
est-neighbor predictor (FIRST-K) achieves 61% accuracy 
(268/440). The accuracies of all predictors are consider- 
ably higher on top scoring subsets of published miRNA. 

[0148] Finally, in order to validate the efficacy and accuracy of 
the DICER-CUT LOCATION DETECTOR 116, a sample of 
novel oligonucleotides detected thereby is preferably se- 
lected, and validated by wet lab. Laboratory results vali- 
dating the efficacy of the DICER-CUT LOCATION DETEC- 



TOR 116 are described hereinbelow with reference to Figs. 
12 -15D. Laboratory confirmation of bioinformatically pre- 
dicted oligonucleotides of the viral GAM group of oligonu- 
cleotides , and several bioinformatically predicted viral 
polynucleotides of the GR group of polynucleotides, is de- 
scribed hereinbelow with reference to Fig. 16. 

[0149] Reference is now made to Fig. 6C which is a simplified 

flowchart illustrating operation of DICER-CUT LOCATION 
DETECTOR 116 (Fig. 2), constructed and operative in ac- 
cordance with a preferred embodiment of the present in- 
vention. The DICER CUT LOCATION DETECTOR 116 
preferably comprises a machine learning computer pro- 
gram module, which is trained to recognize dicer-cut lo- 
cations on known hairpin-shaped miRNA precursors, and 
based on this training, is operable to detect dicer-cut lo- 
cations of novel GAM RNAs (Fig. 1) on GAM FOLDED PRE- 
CURSOR RNAs (Fig. 1). In a preferred embodiment of the 
present invention, the dicer-cut location module prefer- 
ably utilizes machine learning algorithms, such as known 
Support Vector Machine (SVM) and more preferably: 
known Bayesian modeling, Nearest Neighbors, and K- 
nearest-neighbor algorithms. 

[0150] when initially assessing a novel GAM FOLDED PRECURSOR 



RNA, all 19-24 nucleotide long segments thereof are ini- 
tially considered as "potential GAM RNAs", since the dicer- 
cut location is initially unknown. 

[0151] For each such potential GAM RNA, the location of its 5' 
end or the locations of its 5' and 3' ends are scored by at 
least one recognition classifier or predictor. 

[° 152 ] In a preferred embodiment of the present invention, the 
DICER-CUT LOCATION DETECTOR 116 (Fig. 2) may use a 
Support Vector Machine predictor trained on and operat- 
ing on features such as the following: 

[0153] Locations of the 5' and/or 3' ends of the known diced 
miRNAs, which are preferably represented by their re- 
spective distances from the 5 end of the corresponding 
hairpin shaped miRNA precursor. Additionally or alterna- 
tively, the 5' and/or 3' ends of the known diced miRNAs 
are preferably represented by the relationship between 
their locations and the locations of one or more nu- 
cleotides along the hairpin shaped miRNA precursor. Ad- 
ditionally or alternatively, the 5' and/or 3' ends of the 
known diced miRNAs are preferably represented by the 
relationship between their locations and the locations of 
one or more bound nucleotide pairs along the hairpin 
shaped miRNA precursor. Additionally or alternatively, the 



5' and/or 3' ends of the known diced miRNAs are prefer- 
ably represented by the relationship between their loca- 
tions and the locations of one or more mismatched nu- 
cleotide pairs along the hairpin shaped miRNA precursor. 
Additionally or alternatively, the 5' and/or 3' ends of the 
known diced miRNAs are preferably represented by the 
relationship between their locations and the locations of 
one or more unmatched nucleotides along the hairpin 
shaped miRNA precursor. Additionally or alternatively, lo- 
cations of the 5' and/or 3' ends of the known diced miR- 
NAs are preferably represented by their respective dis- 
tances from the loop located at the center of the corre- 
sponding hairpin shaped miRNA precursor; and secondar- 
ily 

[0154] Metrics related to the nucleotide content of the diced 
miRNA and/or of the hairpin shaped miRNA precursor. 

[0155] | n another preferred embodiment of the present inven- 
tion, the DICER-CUT LOCATION DETECTOR 116 (Fig. 2) 
preferably employs an "EDIT DISTANCE" predictor, which 
seeks sequences that are similar to those of known miR- 
NAs, utilizing a Nearest Neighbor algorithm, where a sim- 
ilarity metric between two sequences is a variant of the 
Edit Distance algorithm (Gusfield, 1997). The EDIT DIS- 



TANCE predictor is based on an observation that miRNA 
oligonucleotides tend to form clusters, the members of 
which show marked sequence similarity. 

[0156] | n y e t another preferred embodiment of the present in- 
vention, the DICER-CUT LOCATION DETECTOR 116 (Fig. 2) 
preferably uses a "TWO PHASE" predictor, which predicts 
the dicer-cut location in two distinct phases: (a) selecting 
a double-stranded segment of the GAM FOLDED PRECUR- 
SOR RNA (Fig. 1) comprising the GAM RNA by nave 
Bayesian modeling and (b) detecting which strand of the 
double-stranded segment contains GAM RNA (Fig. 1) by 
employing either nave or by K-nearest-neighbor model- 
ing. K-nearest-neighbor modeling is a variant of the 
'FIRST-K' predictor described hereinbelow, with parame- 
ters optimized for this specific task. The 'TWO PHASE' pre- 
dictor may be operated in two modes: either utilizing only 
the first phase and thereby producing two alternative 
dicer-cut location predictions, or utilizing both phases 
and thereby producing only one final dicer-cut location. 

[° 157 ] In still another preferred embodiment of the present in- 
vention, the DICER-CUT LOCATION DETECTOR 116 
preferably uses a "FIRST-K" predictor, which utilizes a K- 
nearest-neighbor algorithm. The similarity metric between 



any two sequences is 1- E/L , where Lisa parameter, 
preferably 8-10 and E is the edit distance between the two 
sequences, taking into account only the first L nucleotides 
of each sequence. If the K-nearest-neighbor scores of two 
or more locations on the GAM FOLDED PRECURSOR RNA 
(Fig. 1) are not significantly different, these locations are 
further ranked by a Bayesian model, similar to the one de- 
scribed hereinabove. 

[0158] T he TWO PHASE and FIRST-K predictors preferably are 

trained on and operate on features such as the following: 

[0159] Locations of the 5' and/or 3' ends of the known diced 
miRNAs, which are preferably represented by their re- 
spective distances from the 5' end of the corresponding 
hairpin shaped miRNA precursor. Additionally or alterna- 
tively, the 5' and/or 3' ends of the known diced miRNAs 
are preferably represented by the relationship between 
their locations and the locations of one or more nu- 
cleotides along the hairpin shaped miRNA precursor. Ad- 
ditionally or alternatively, the 5' and/or 3' ends of the 
known diced miRNAs are preferably represented by the 
relationship between their locations and the locations of 
one or more bound nucleotide pairs along the hairpin 
shaped miRNA precursor. Additionally or alternatively, the 



5' and/or 3' ends of the known diced miRNAs are prefer- 
ably represented by the relationship between their loca- 
tions and the locations of one or more mismatched nu- 
cleotide pairs along the hairpin shaped miRNA precursor. 
Additionally or alternatively, the 5' and/or 3' ends of the 
known diced miRNAs are preferably represented by the 
relationship between their locations and the locations of 
one or more unmatched nucleotides along the hairpin 
shaped miRNA precursor. Additionally or alternatively, lo- 
cations of the 5' and/or 3' ends of the known diced miR- 
NAs are preferably represented by their respective dis- 
tances from the loop located at the center of the corre- 
sponding hairpin shaped miRNA precursor; and secondar- 
ily Metrics related to the nucleotide content of the diced 
miRNA and/or of the hairpin shaped miRNA precursor. 
[0160] | n accordance with an embodiment of the present inven- 
tion scores of two or more of the abovementioned classi- 
fiers or predictors are integrated, yielding an integrated 
score for each "potential GAM RNA". As an example, Fig. 
6C illustrates integration of scores from two classifiers, a 
3' end recognition classifier and a 5' end recognition clas- 
sifier, the scores of which are integrated to yield an inte- 
grated score. Most preferably, the INTEGRATED SCORE of 



Fig. 6C preferably implements a "best-of-breed" approach 
employing a pair of classifiers and accepting only "poten- 
tial GAM RNAs" that score highly on one of the above 
mentioned "EDIT DISTANCE", or "TWO-PHASE" predictors. 
In this context, "high scores" means scores which have 
been demonstrated to have low false positive value when 
scoring known miRNA oligonucleotides. Alternatively, the 
INTEGRATED SCORE may be derived from operation of 
more or less than two classifiers. 

[° 161 ] The INTEGRATED SCORE is evaluated as follows: (a) the 
"potential GAM RNA" having the highest score is prefer- 
ably taken to be the most probable GAM RNA, and (b) if 
the integrated score of this most probable GAM RNA is 
higher than a pre-defined threshold, then the most prob- 
able GAM RNA is accepted as a PREDICTED GAM RNA. 
Preferably, this evaluation technique is not limited to the 
highest scoring potential GAM RNA. 

[0162] as a final optional stage we may filter PREDICTED GAM 

RNAs with low complexity, which have high probability to 
be part of repeated element in the DNA (for example 
ATATATA sequences). For each PREDICTED GAM RNA se- 
quence we count the number of occurrences of each 2 
nucleotides combination (AA, AT, AC), we filter out se- 



quences where the sum of the most 2 probable combina- 
tion is higher then 9. Using this criteria we filter 2% of the 
known miRNA and around 30% of the PREDICTED GAM 
RNAs. 

[0163] Reference is now made to Fig. 7A which is a simplified 

block diagram of a preferred implementation of the TAR- 
GET GENE BINDING-SITE DETECTOR 118 described herein- 
above with reference to Fig. 2. The goal of the TARGET 
GENE BINDING-SITE DETECTOR 118 is to detect one or 
more binding sites such as BINDING SITE I, BINDING SITE II 
and BINDING SITE III (Fig. 1) located in 3' untranslated re- 
gions of the mRNA of a known gene, the nucleotide se- 
quence of which binding sites is partially or fully comple- 
mentary to a GAM RNA, thereby determining that the 
above mentioned known gene is a target gene thereof. 

[° 1 64] The TARGET GENE BINDING-SITE DETECTOR 118 (Fig. 2) 

receives a plurality of DICER-CUT SEQUENCES FROM HAIR- 
PIN STRUCTURES 140 (Fig. 6A), and a plurality of POTEN- 
TIAL TARGET GENE SEQUENCES 142 which are derived 
from SEQUENCED DNA DATA 104 (Fig. 2). 

[0165] TARGET GENE BINDING-SITE DETECTOR TRAINING & VALI- 
DATION FUNCTIONALITY 128 (Fig. 3) is operative to train 
the TARGET GENE BINDING-SITE DETECTOR 118 on known 



miRNAs and their respective target genes and to build a 
background model for evaluation of a P value for the TAR- 
GET GENE BINDING-SITE DETECTOR118 resultslt Construct 
the model by analyzing both heuristically and computa- 
tionally the results of the TARGET GENE BINDING-SITE DE- 
TECTOR 118. 

[0 1 66] Following operation of TARGET GENE BINDING-SITE DE- 
TECTOR TRAINING & VALIDATION FUNCTIONALITY 128 
(Fig. 3), the TARGET GENE BINDING-SITE DETECTOR 118 is 
operative to detect a plurality of POTENTIAL NOVEL TAR- 
GET GENES HAVING BINDING-SITE/S 144 the nucleotide 
sequence of which is partially or fully complementary to 
that of each of the plurality of DICER-CUT SEQUENCES 
FROM HAIRPIN STRUCTURES 140. Preferred operation of 
the TARGET GENE BINDING-SITE DETECTOR 118 is further 
described hereinbelow with reference to Fig. 7B. 

[0167] Reference is now made to Fig. 7B which is a simplified 

flowchart illustrating a preferred operation of the TARGET 
GENE BINDING-SITE DETECTOR 118 of Fig. 2. In a pre- 
ferred embodiment of the present invention, the TARGET 
GENE BINDING-SITE DETECTOR 118 employs a match 
search tool .described hereinbelow, in order to compare 
the nucleotide sequence of each of the plurality of DICER- 



CUT SEQUENCES FROM HAIRPIN STRUCTURES 140 (Fig. 
6A), to the POTENTIAL TARGET GENE SEQUENCES 142 (Fig. 
7 A), such as 3 untranslated regions of known mRNAs, in 
order to find crude potential matches. Alternatively, a se- 
quence comparison algorithm such as BLAST can be used. 

[0168] The match search tool compare the DICER-CUT SE- 
QUENCES FROM HAIRPIN STRUCTURES 140 (Fig. 7A) to PO- 
TENTIAL TARGET GENE SEQUENCES 142 (Fig. 7A). the 
methodology of the match search tool is to match Nu- 
cleotides 1:8 of DICER-CUT SEQUENCES FROM HAIRPIN 
STRUCTURES 140 (Fig. 7A) to POTENTIAL TARGET GENE 
SEQUENCES 142 (Fig. 7A).The results are preferably fil- 
tered to allow at most 1 edit distance between the DICER- 
CUT SEQUENCES FROM HAIRPIN STRUCTURES 140 (Fig. 7A) 
complementary sequence and POTENTIAL TARGET GENE 
SEQUENCES 142 (Fig. 7A). 

[0169] | n a preferred implementation of the match search tool, it 
first calculates two preprocess data structures: the first, 
store for every 9nts sequences all its appearances in the 
POTENTIAL TARGET GENE SEQUENCES 142 (Fig. 7A). 

[0170] The second is a map of all the 8nts sequence to 9nts se- 
quences that can be accepted by adding one nucleotide 
and performing one or less editing actions to the 8nts 



complementary sequence. This two data structures help to 
improve performance and to perform the binding site fil- 
tering stage as an integral part of the search. The match 
search tool in this implementation returns for every 
DICER-CUT SEQUENCES FROM HAIRPIN STRUCTURES 140 
(Fig. 7 A all the appearances of the 9nts sequence in the 
POTENTIAL TARGET GENE SEQUENCES 142 (Fig. 7A) that 
are calculated using the second, describe hereinabove, 
data structure from its first 8nts. The 9nts sequences ap- 
pearances are the potentially 3 of the POTENTIAL TARGET 
GENE SEQUENCES 142 (Fig. 7A) binding site of each 
DICER-CUT SEQUENCES FROM HAIRPIN STRUCTURES 140 
(Fig. 7A 

[0171] Next the binding sites are expanded, and determinations 
are made whether if nucleotide sequences immediately 
upstream to the binding sites 3 part matchs with the 3' 
end of the DICER-CUT SEQUENCES FROM HAIRPIN STRUC- 
TURES. The alignment of the binding site 5 part is made 
by the alignment algorithm with specific weighting pa- 
rameters describes hereinbelow with reference to Fig. 5). 
Alternatively the alignment can include also the 3 part of 
the binding site. 

[0172] Th e alignment algorithm (implemented here) with specific 



weighting parameters. The weighting parameters are 
based on analyze of known miRNA binding sites 
(references: Wightman et al. Cell 1993, Moss Cell 1997, 
Reinhart Nature 2000, Abrahante Developmental Cell 
2003 and Lin Developmental Cell 2003, Brennecke Cell 
2003, Stark PLOS 2003, lewis et al. Cell 2003). It was 
heuristically found to best distinguish between known 
miRNA binding sites and sequences which are known not 
to be miRNA binding sites. 
[0173] According to a preferred embodiment of the present in- 
vention a good matching at 5' end of the miRNA 
(nucleotides 1-8), and at 3' end of the miRNA, (the last 9 
nucleotides) is required. The matching at 5' end is crucial 
for the binding while the 3' end binding is less important 
but can compensate on few mismatches at the 5'. The nu- 
cleotides at the middle of the miRNA have no punishment 
for any editing operation (reference: Doench et al. Gene 
and Development 2004). In this preferred embodiment the 
alignment of the 5 end of the miRNA (3 of the binding 
site) is a result of the match search tool, while the middle 
and 3 miRNA (5 of the binding site) is done by the align- 
ment algorithm. The alignment is a combination of both. 
Alternatively, the alignment is done for all the miRNA and 



all the binding site sequence. The alignment algorithm 
finds the lowest free-energy spatial structure and heuris- 
tic constrain that were found to distinct the known bind- 
ing sites in the best way. The alignment yields 2 scores 
for 5' end and 3' end. Each score represents the inverse of 
its probability to obtain such or higher score by running 
random sequences, with 22 nucleotides length and same 
base composition as known miRNA, with randomly chosen 
3 un translated region sequences. The final score is the 
multiplicity of the two scores. 
[0174] The binding site detector, then analyze the spatial struc- 
ture of the miRNA with all its potential target gene 3 un 
translated region. It assesses the list of binding sites of 
each miRNA - Target gene pair with their positions and 
respective scores. The TARGET GENE BINDING-SITE DE- 
TECTOR 118 (Fig. 7A) then uses the fact that Many of the 
known binding sites are clustered, to evaluate the P value 
of obtaining cluster of few binding sites on the same tar- 
get gene 3 un translated region in the following way: It 
scans different scores threshold, and calculate for each 
threshold the number and positions of possible binding 
sites with score above the threshold. It then gets from a 
preprocess calculated background matrixes, describe 



hereinbelow, a P value for each threshold and number and 
positions of binding sites combination. The output score 
for each miRNA and target gene pair is the minimal P 
value, normalized with the number of thresholds trying 
using Bernoulli distribution. A preference of low P value 
pairs is made. 

[0175] As mentioned hereinabove, for each target gene a prepro- 
cess calculated background matrixes were build. The ma- 
trixes includes rows for each number of miRNA binding 
site (In the preferred embodiment the matrix include 7 
rows for 0 to 6 binding site.), and columns for each dif- 
ferent score threshold(ln the preferred embodiment the 
matrix include 5 columns for 5 different threshold). Each 
matrix cell, corresponding to specific number of binding 
site and threshold, was set to be the probability to get 
equal or more binding sites with equal or more score by 
random sequences with 22 nucleotides length with the 
same letter distribution as known miRNAs (29.5% T, 24.5% 
A, 25% G and 21% C). Those probabilities are calculated by 
running the above procedure for 10000 random se- 
quences. The P value can be estimated as the number of 
random sequences obeys the matrix cell requirement di- 
vided by the total number of random sequences (10000). 



In the preferred embodiment 2 matrixes are calculated. 
The P values are of the second matrix are calculated under 
constrain that at least two of the binding site position are 
under a heuristically determine constant value. The sec- 
ond matrix values are calculated without this constrain. 
The TARGET GENE BINDING-SITE DETECTOR uses the sec- 
ond matrix if the binding sites position agree with the 
constrain, else it uses the first. In an alternative embodi- 
ment only one matrix is calculated without any constrain 
on the binding sites positions. 

[° 176 ] Test preformed using the TARGET GENE BINDING-SITE 
DETECTOR showed that all of the known miRNA target 
genes are found using this algorithm with a P value of less 
than 0.5 %. Running known miRNA against 3400 potential 
3' UTR of target genes sequences yield on average 32 tar- 
get genes for each miRNA with P value less than 0.5 %, 
while background sequences, as well as inverse or com- 
plement sequence of known miRNA (which preserve their 
high order sequence statistics) found, as excepted, on av- 
erage 17 target genes. This reflects the algorithm perfor- 
mance to detect real target genes with 47% accuracy. 

[0177] Scoring of all potential binding sites are calculated and P 
value for obtaining such binding sites list for DICER-CUT 



SEQUENCE FROM HAIRPIN STRUCTURE and POTENTIAL 
TARGET GENE SEQUENCE pair is evaluate using the back- 
ground matrix. In accordance with another preferred em- 
bodiment of the present invention, binding sites are 
searched by a reversed process. Sequences of K 
(preferably 22) nucleotides of a untranslated regions of a 
target gene are assessed as potential binding sites. A se- 
quence comparison algorithm, such as BLAST or EDIT DIS- 
TANCE varinet, is then used to search elsewhere in the 
genome for partially or fully complementary sequences 
which are found in known miRNA oligonucleotides or 
computationally predicted GAM oligonucleotides. Only 
complementary sequences, which meet predetermined 
spatial structure and free energy criteria as described 
hereinabove, are accepted. Clustered binding sites are 
strongly preferred and potential binding sites and poten- 
tial GAM oligonucleotides which occur in evolutionarily 
conserved genomic sequences are also preferred. Scoring 
of candidate binding sites takes into account free energy 
and spatial structure of the binding site complexes, as 
well as the aforesaid preferences. 
[0178] UTR of GAM viral target genes were preferably extracted 
directly from annotation of UTR records. Alternatively, 



UTR of GAM viral target genes were preferably extracted 
by taking the sequences spanned from last coding posi- 
tion to the 3 end of the mRNA sequence annotation. Al- 
ternatively, UTR of GAM viral target genes were preferably 
extracted by taking 400 nucleotides downstream to the 
endcoding region of the mRNA sequence. All of above- 
mention methods were applied on complete viral genomes 
data in GeneBank format from the NCBI RefSeq database, 
version 18-Jan-2004 

(ftp://ftp.ncbi.nih.gov/refseq/release/viral). 
[0179] Reference is now made to Fig. 8 which is a simplified 

flowchart illustrating a preferred operation of the FUNC- 
TION & UTILITY ANALYZER 120 described hereinabove 
with reference to Fig. 2. The goal of the FUNCTION & 
UTILITY ANALYZER 120 is to determine if a potential tar- 
get gene is in fact a valid clinically useful target gene. 
Since a potential novel GAM oligonucleotide binding a 
binding site in the UTR of a target gene is understood to 
inhibit expression of that target gene, and if that target 
gene is shown to have a valid clinical utility, then in such a 
case it follows that the potential novel oligonucleotide it- 
self also has a valid useful function which is the opposite 
of that of the target gene. 



[0180] T he FUNCTION & UTILITY ANALYZER 120 preferably re- 
ceives as input a plurality of POTENTIAL NOVEL TARGET 
GENES HAVING BINDING-SITE/S 144 (Fig. 7A), generated 
by the TARGET GENE BINDING-SITE DETECTOR 118 (Fig. 
2). Each potential oligonucleotide is evaluated as follows: 
First, the system checks to see if the function of the po- 
tential target gene is scientifically well established. Prefer- 
ably, this can be achieved bioinformatically by searching 
various published data sources presenting information on 
known function of proteins. Many such data sources exist 
and are published as is well known in the art. Next, for 
those target genes the function of which is scientifically 
known and is well documented, the system then checks if 
scientific research data exists which links them to known 
diseases. For example, a preferred embodiment of the 
present invention utilizes the OMIM(TM) (Hamosh et al, 
2002) database published by NCBI, which summarizes re- 
search publications relating to genes which have been 
shown to be associated with diseases. Finally, the specific 
possible utility of the target gene is evaluated. While this 
process too may be facilitated by bioinformatic means, it 
might require manual evaluation of published scientific 
research regarding the target gene, in order to determine 



the utility of the target gene to the diagnosis and or treat- 
ment of specific disease. Only potential novel oligonu- 
cleotides, the target genes of which have passed all three 
examinations, are accepted as novel oligonucleotide. 

[0181] Reference is now made to Fig. 9, which is a simplified dia- 
gram describing each of a plurality of novel bioinformati- 
cally detected regulatory polynucleotide referred to here 
as Genomic Record (GR) polynucleotide, which encodes an 
x operon-like x cluster of novel microRNA-like oligonu- 
cleotides, each of which in turn modulates expression of 
at least one target gene; the function and utility of which 
at least one target gene is known in the art. 

[0182] GR PRECURSOR is a novel bioinformatically detected regu- 
latory, non protein coding polynucleotide. The method by 
which GR PRECURSOR is detected is described hereinabove 
with additional reference to Figs. 1-9. 

[0183] gr PRECURSOR is preferably encoded by a viral genome 
and contains cluster of novel viral oligonucleotides, which 
preferably binds to a human target genes or to the virus 
self-genes. Alternatively or additionally, GR PRECURSOR is 
encoded by the human genome and contains cluster of 
novel human oligonucleotides, which preferably binds to 
viral target genes or to human self-genes. 



[0 1 84] GR PRECURSOR encodes GR PRECURSOR RNA, an RNA 
molecule, typically several hundred to several thousand 
nucleotides long. GR PRECURSOR RNA folds spatially, 
forming GR FOLDED PRECURSOR RNA. It is appreciated 
that GR FOLDED PRECURSOR RNA comprises a plurality of 
what is known in the art as x hairpin x structures. These 
x hairpin x structures are due to the fact that the nu- 
cleotide sequence of GR PRECURSOR RNA comprises a 
plurality of segments, the first half of each such segment 
having a nucleotide sequence which is at least a partial or 
accurate reverse-complement sequence of the second half 
thereof, as is well known in the art. 

[° 1 85] GR FOLDED PRECURSOR RNA is naturally processed by 

cellular enzymatic activity into a plurality of separate GAM 
precursor RNAs, herein schematically represented by 
GAM1 FOLDED PRECURSOR RNA through GAM 3 FOLDED 
PRECURSOR RNA, each of which GAM folded precursor 
RNAs being a hairpin shaped RNA segment, correspond- 
ing to GAM FOLDED PRECURSOR RNA of Fig.l. 

[0186] The above mentioned GAM folded precursor RNAs are 
diced by DICER COMPLEX of Fig.l, yielding respective 
short RNA segments of about 22 nucleotides in length, 
schematically represented by GAM1 RNA through GAM3 



RNA, each of which GAM RNAs corresponding to GAM RNA 
of Fig.l.GAMl RNA, GAM 2 RNA and GAM 3 RNA, each bind 
complementarily to binding sites located in untranslated 
regions of respective target genes, designated GAM1 
TARGET RNA, GAM2 TARGET RNA and GAM 3 TARGET RNA, 
respectively, which target binding site corresponds to a 
target binding site such as BINDING SITE I, BINDING SITE II 
or BINDING SITE III of Fig.l. This binding inhibits transla- 
tion of the respective target proteins designated GAM1 
TARGET PROTEIN, GAM 2 TARGET PROTEIN and GAM 3 
TARGET PROTEIN respectively. 

[0187] it i S appreciated that specific functions, and accordingly 
utilities of GR polynucleotide are correlated with and may 
be deduced from the identity of the target genes, which 
are inhibited by GAM RNAs comprised in the 
x operon-like^ cluster of GR polynucleotide, schematically 
represented by GAM1 TARGET PROTEIN through GAM 3 
TARGET PROTEIN. The function of these target genes is 
elaborated in Table 8, hereby incorporated herein. 

[0188] Reference is now made to Fig. 10 which is a block diagram 
illustrating different utilities of oligonucleotide of the 
novel group of oligonucleotides of the present invention 
referred to here as GAM oligonucleotides and GR polynu- 



cleotides. The present invention discloses a first plurality 
of novel oligonucleotides referred to here as GAM 
oligonucleotides and a second plurality of operon-like 
polynucleotides referred to here as GR polynucleotides, 
each of the GR polynucleotide encoding a plurality of GAM 
oligonucleotides. The present invention further discloses a 
very large number of known target genes, which are 
bound by, and the expression of which is modulated by 
each of the novel oligonucleotides of the present inven- 
tion. Published scientific data referenced by the present 
invention provides specific, substantial, and credible evi- 
dence that the above mentioned target genes modulated 
by novel oligonucleotides of the present invention, are as- 
sociated with various diseases. Specific novel oligonu- 
cleotides of the present invention, target genes thereof 
and diseases associated therewith, are described herein- 
below with reference to Tables 1 through 10 . It is there- 
fore appreciated that a function of GAM oligonucleotides 
and GR polynucleotides of the present invention is modu- 
lation of expression of target genes related to known viral 
diseases, and that therefore utilities of novel oligonu- 
cleotides of the present invention include diagnosis and 
treatment of the above mentioned diseases. 



[0189] pig. 10 describes various types of diagnostic and thera- 
peutic utilities of novel oligonucleotides of the present in- 
vention. A utility of novel oligonucleotide of the present 
invention is detection of GAM oligonucleotides and of GR 
polynucleotides. It is appreciated that since GAM oligonu- 
cleotides and GR polynucleotides modulate expression of 
disease related target genes, that detection of expression 
of GAM oligonucleotides in clinical scenarios associated 
with said viral diseases is a specific, substantial and credi- 
ble utility. Diagnosis of novel oligonucleotides of the 
present invention may preferably be implemented by RNA 
expression detection techniques, including but not limited 
to biochips, as is well known in the art. Diagnosis of ex- 
pression of oligonucleotides of the present invention may 
be useful for research purposes, in order to further un- 
derstand the connection between the novel oligonu- 
cleotides of the present invention and the above men- 
tioned related viral diseases, for disease diagnosis and 
prevention purposes, and for monitoring disease 
progress. 

[0190] Another utility of novel oligonucleotides of the present in- 
vention is anti-GAM therapy, a mode of therapy which al- 
lows up regulation of a viral disease-related target gene 



of a novel GAM oligonucleotide of the present invention, 
by lowering levels of the novel GAM oligonucleotide which 
naturally inhibits expression of that target gene. This 
mode of therapy is particularly useful with respect to tar- 
get genes which have been shown to be under-expressed 
in association with a specific viral disease. Anti-GAM ther- 
apy is further discussed hereinbelow with reference to 
Figs. 11A and 11B. 
[0191] a further utility of novel oligonucleotides of the present 
invention is GAM replacement therapy, a mode of therapy 
which achieves down regulation of a viral disease related 
target gene of a novel GAM oligonucleotide of the present 
invention, by raising levels of the GAM which naturally in- 
hibits expression of that target gene. This mode of ther- 
apy is particularly useful with respect to target genes 
which have been shown to be over-expressed in associa- 
tion with a specific viral disease. GAM replacement ther- 
apy involves introduction of supplementary GAM products 
into a cell, or stimulation of a cell to produce excess GAM 
products. GAM replacement therapy may preferably be 
achieved by transfecting cells with an artificial DNA 
molecule encoding a GAM which causes the cells to pro- 
duce the GAM product, as is well known in the art. 



[0192] yet a further utility of novel oligonucleotides of the 

present invention is modified CAM therapy. Disease con- 
ditions are likely to exist, in which a mutation in a binding 
site of a GAM RNA prevents natural GAM RNA to effec- 
tively bind inhibit a viral disease related target gene, 
causing up regulation of that target gene, and thereby 
contributing to the disease pathology. In such conditions, 
a modified GAM oligonucleotides is designed which effec- 
tively binds the mutated GAM binding site, i.e. is an effec- 
tive anti-sense of the mutated GAM binding site, and is 
introduced in disease effected cells. Modified GAM ther- 
apy is preferably achieved by transfecting cells with an ar- 
tificial DNA molecule encoding the modified GAM which 
causes the cells to produce the modified GAM product, as 
is well known in the art. 

[0193] Reference is now made to Figs. 11A and 11B, simplified 
diagrams which when taken together illustrate anti-GAM 
therapy mentioned hereinabove with reference to Fig. 10. 
A utility of novel GAMs of the present invention is anti- 
GAM therapy, a mode of therapy which allows up regula- 
tion of a viral disease-related target gene of a novel GAM 
of the present invention, by lowering levels of the novel 
GAM which naturally inhibits expression of that target 



gene. Fig. 11A shows a normal GAM inhibiting translation 
of a target gene by binding of CAM RNA to a BINDING SITE 
found in an untranslated region of GAM TARGET RNA, as 
described hereinabove with reference to Fig. 1. 
[0194] pig. 11B shows an example of anti-GAM therapy. ANTI- 
GAM RNA is short artificial RNA molecule the sequence of 
which is an anti-sense of GAM RNA. Anti-GAM treatment 
comprises transfecting diseased cells with ANTI-GAM 
RNA, or with a DNA encoding thereof. The ANTI-GAM RNA 
binds the natural GAM RNA, thereby preventing binding of 
natural GAM RNA to its BINDING SITE. This prevents natu- 
ral translation inhibition of GAM TARGET RNA by GAM 
RNA, thereby up regulating expression of GAM TARGET 
PROTEIN. 

[0195] it is appreciated that anti-GAM therapy is particularly use- 
ful with respect to target genes which have been shown to 
be under-expressed in association with a specific viral 
disease. 

[0196] Furthermore, anti-GAM therapy is particularly useful, 

since it may be used in situations in which technologies 
known in the art as RNAi and siRNA can not be utilized. As 
in known in the art, RNAi and siRNA are technologies 
which offer means for artificially inhibiting expression of a 



target protein, by artificially designed short RNA segments 
which bind complementarily to mRNA of said target pro- 
tein. However, RNAi and siRNA can not be used to directly 
up regulate translation of target proteins. 
[0197] Reference is now made to Fig. 12A, which is a histogram 
representing the distribution of known miRNA oligonu- 
cleotides and that of hairpin structures extracted from ex- 
pressed human genome sequences with respect to their 
HAIRPIN DETECTOR score (this analysis preferably uses 
the previous version of the Hairpin Detector program). 
The known miRNA oligonucleotides set is taken from 
RFAM database, Release2.1 and includes 440 miRNA 
oligonucleotides from H.sapienas, M.musculus, C.elegans, 
C.brigassae and D.melanogaster. Folding of expressed 
genome sequences taken from public databases of ESTs 
(Unigene-NCBI and TIGR) identified 342,882 hairpin 
structures. -154,000 out of the 342,882 hairpin structures 
did not pass the filter of being identified as hairpins in 
several secondary structure folding versions of the given 
genomic sequence , as described hereinabove with refer- 
ence to Fig. 5B, and hence did not receive a Hairpin detec- 
tor score. Furthermore, ~133,000 hairpin structures did 
not pass the filter of minimum score of the DICER-CUT 



LOCATION DETECTOR 116 (Fig. 2) (those -287,000 hair- 
pin structures are not represented in the histogram). Hair- 
pin structures are considered as miRNA -like precursor 
oligonucleotides here referred to as GAM precursor, if 
their Hairpin detector score is above 0.3. Thus, the GAM 
precursor set is comprised of ~40,000 hairpin structures, 
of those ~5100 received a high Hairpin detector score 
(>=0.7). These are much higher numbers than those of 
the known miRNA oligonucleotides and of the upper 
bound of ~255 human miRNA oligonucleotides, estimated 
by Bartel et al. (Science, 299, 1540, March 2003). Of the 
reference set that pass the above filter (408/440), 284 
(69%) received a high Hairpin detector score (>=0.7). 
[0198] Reference is now made to Fig. 12B, which is a table sum- 
marizing laboratory validation results that validate efficacy 
of the BIOINFORMATIC OLIGONUCLEOTIDE DETECTION 
ENGINE 100 ( Fig. 2). In order to assess efficacy of the 
BIOINFORMATIC OLIGONUCLEOTIDE DETECTION ENGINE 
100, novel oligonucleotides predicted thereby are prefer- 
ably filtered using the scores of HAIRPIN DETECTOR 114 
(Fig. 2) and DICER-CUT LOCATION DETECTOR 116 (Fig. 2) 
as follows: covering 85% (second column) of published 
hairpin our algorithm enable to filter 95% of background 



(third column), leaving only dozens of predicted hairpins 
with relatively high probability of being a real hairpin car- 
rying a 'diced' miRNA for lab validation. 

[0199] Group A (selected for the present invention, first column): 
The score of the HAIRPIN-DETECTOR is above 0, the over- 
all score of the two-phased predictor is above 0.55, and 
the score of the second phase of the two-phased predic- 
tor is above 0.75, or the score of the EDIT-DISTANCE pre- 
dictor is equal or above 17. In this group, one Dicer cut 
location is predicted for each hairpin. 

[0200] Sample of novel bioinformatically predicted human GAMs 
are sent to the laboratory for validation (fourth column), 
and the number (fifth column) and percent (sixth column) 
of successful validation of predicted human CAM is noted 
for each of the groups, as well as overall (bottom line). 
The information about the above lab validation is based 
on previous patents group selections, there is a high cor- 
relation between the old group A predicted hairpin and 
the present one indicate that the update results should be 
approximately the same or better, due to improvement in 
the predictions. 

[0201] it is appreciated that the present invention comprises 

15484 novel GAM oligonucleotides, which fall into group 



A, and that the BIOINFORMATIC OLIGONUCLEOTIDE DE- 
TECTION ENGINE 100 ( Fig. 2) is substantiated by a group 
of 37novel human GAM oligonucleotides validated by lab- 
oratory means, out of 101 human GAM oligonucleotides 
which were tested in the lab, resulting in validation of an 
overall 31% accuracy. The selected group demonstrated 
37% accuracy. Pictures of test-results of specific human 
GAM oligonucleotides as well as the methodology used for 
validating the expression of predicted oligonucleotides 
are elaborated hereinbelow with reference to Fig. 13. 

[0202] | t j S further appreciated that failure to detect a predicted 
oligonucleotide in the lab does not necessarily indicate a 
mistaken bioinformatic prediction. Rather, it may be due 
to technical sensitivity limitation of the lab test, or be- 
cause the predicted oligonucleotides are not expressed in 
the tissue examined, or at the development phase tested. 

[0203] it is still further appreciated that in general these findings 
are in agreement with the expected bioinformatic accu- 
racy, as describe hereinabove with reference to Fig. 6B: 
assuming 80% accuracy of the HAIRPIN DETECTOR 114 
and 80% accuracy of the DICER-CUT LOCATION DETECTOR 
116 and 80% accuracy of the lab validation, this would re- 
sult in 50% overall accuracy of the predicted oligonu- 



cleotide validated in the lab. 
[0204] Reference is now made to Fig. 13A which is a picture of 
laboratory results validating the expression of 43 novel 
human GAM oligonucleotides detected by the BIOINFOR- 
MATIC OLIGONUCLEOTIDE DETECTION ENGINE 100 (Fig. 
2). 

[0205] Reference is now made to Fig. 13A and Fig. 13B which are 
pictures and a summary table of laboratory results vali- 
dating the expression of 43 novel human GAM oligonu- 
cleotides detected by the BIOINFORMATIC OLIGONU- 
CLEOTIDE DETECTION ENGINE 100. In each row in Fig. 
13A, pictures of several oligonucleotides validated by hy- 
bridization of Polymerase Chain Reaction (PCR)-product 
southern-blots, are provided, each corresponding to a 
specific GAM oligonucleotide, as elaborated hereinbelow. 
To test our validation method, we used a reference set of 
8 known human miRNA oligonucleotides, as a blind test to 
our laboratory. These PCR-product hybridization pictures 
are designated 1 through 8 for the reference set known 
miRNA oligonucleotides; and 9 through 51 for predicted 
GAM oligonucleotides. 

[0206] | n each PCR hybridization picture, 2 lanes are seen: the 

test lane, designated "+" and the control lane, designated 



The control reaction contained all the components of 
the test reaction except library template. For convenience 
of viewing the results, all PCR-product hybridization pic- 
tures of Fig. 22A have been shrunk x4 vertically. It is ap- 
preciated that for each of the tested GAM oligonucleotides 
a clear hybridization band appears in the test ("+") lane, 
but not in the control ("-") lane. 

[0207] Specifically, Fig. 13A shows pictures of PCR-product hy- 
bridization validation by southern-blot, the methodology 
of which is described hereinbelow, to the following novel 
human GAM oligonucleotides (RosettaGenomics Ltd. 
Nomenclature, 'A 1 and 'B' referred to the Dicer Cut Loca- 
tion as described hereinabelow with reference to the de- 
scription of large tables: 

[0208] (l)hsa-MIR-21;(2)hsa-MIR-27b; (3)hsa-MIR-186; 
(4)hsa-MIR-93; (5)hsa-MIR-26a ; (6)hsa-MIR-191; 
(7)hsa-MIR-31; (8)hsa-MIR-92; (9) GAM3418-A (later 
published by other researchers as hsa-MIR23); (10) 
GAM4426-A; (11) GAM281-A; (12) GAM7553-A; (13) 
GAM5385-A; (14) GAM2608-A; (15) GAM1032-A; (16) 
GAM3431-A; (17) GAM7933-A; (18) GAM3298-A.; (19) 
GAM7080-A; (20)GAM895-A.; (21) GAM3770.1; (22) 
GAM337162-A; (23) GAM 8678-A; (24) GAM2033-A; 



(25)GAM7776-A; (26) GAM8145-A; (27) GAM25-A; (28) 
GAM7352.1; (29) GAM337624-A; (30) GAM1479-A; (31) 
GAM2270-A; (32) GAM7591-A; (33) GAM8285-A; (34) 
GAM6773-A; (35) GAM336818-A; (36) GAM336487-A; 
(37) GAM337620-A; (38) GAM336809-A; (39) 
GAM5346-A; (40) GAM8554-A; (41)GAM2071-A; (42) 
GAM7957-A; (43) GAM391-A; (44) GAM6633-A; (45) 
GAM19; (46) GAM8358-A; (47) GAM3229-A; an)GAM 
7052-A; (49) GAM3027-A; (50) GAM21and (51) GAM 
oligonucleotide similar to mmu-MIR-30e. 
[0209] The next validated GAM oligonucleotides are highly simi- 
lar or identical to known mouse-miRNA oligonucleotides: 
GAM3027-A, similar to mmu-MIR-29c; GAM21, similar to 
mmu-MIR-130b; and GAM oligonucleotide which is highly 
similar to mmu-MIR-30e (picture number 51). In addition 
to the PCR- product hybridization detection, the following 
GAMs were also cloned and sequenced : GAM3418-A , 
GAM5385-A, GAM1032-A, GAM3298-A, GAM7080-A, 
GAM1338-A, GAM7776-A, GAM25-A, GAM337624-A, 
GAM1479-A, GAM6773-A, GAM336818-A, 
GAM336487-A, GAM337620-A, GAM336809-A, 
GAM3027-A, GAM21, and GAM oligonucleotide similar to 
mmu-MIR-30e (picture number 51). Furthermore, the fol- 



lowing GAM oligonucleotides were sequenced directly 
from the ligation reaction by the method described here- 
inbelow under LIGATION-PCR DIAGNOSTIC METHOD: 
GAM4426-A, GAM7553-A, GAM2270-A, and 
GAM7591-A. 

[0210] | n order to validate the expression of predicted novel 
GAMs and assuming that these novel GAM oligonu- 
cleotides are probably expressed at low concentrations, a 
PCR product cloning approach was set up through the fol- 
lowing strategy: two types of cDNA libraries designated 
"One tailed" and "Ligation" were prepared from frozen 
HeLa S100 extract (4c Biotech, Belgium) size fractionated 
RNA. Essentially, Total S100 RNA was prepared through an 
SDS-Proteinase K incubation followed by an acid Phenol- 
Chloroform purification and Isopropanol precipitation. Al- 
ternatively, total HeLa RNA was also used as starting ma- 
terial for these libraries. 

[021 1] Fractionation was done by loading up to 500g per YM100 
Amicon Microcon column (Millipore) followed by a 500g 
centrifugation for 40 minutes at 4C. Flow through 
"YM100"RNA consisting of about of the total RNA was 
used for library preparation or fractionated further by 
loading onto a YM30 Amicon Microcon column (Millipore) 



followed by a 13,500g centrifugation for 25 minutes at 
4C. Flowthrough "YM30" was used for library preparation 
as is and consists of less than 0.5% of total RNA. For the 
both the "ligation" and the "One-tailed" libraries, RNA was 
dephosphorilated and ligated to an RNA (lowercase)-DNA 
(UPPERCASE) hybrid 5-phosphorilated, 3 idT blocked 
3-adapter (5-P-uuuAACCGCATCCTTCTC-idT-3 Dharma- 
con # P-002045-01-05) (as elaborated in Elbashir et al., 
Genes Dev. 15:188-200 (2001)) resulting in ligation only 
of RNase III type cleavage products. 3-Ligated RNA was 
excised and purified from a half 6%, half 13% polyacry- 
lamide gel to remove excess adapter with a Nanosep 0.2M 
centrifugal device (Pall) according to instructions, and 
precipitated with glycogen and 3 volumes of Ethanol. Pel- 
let was resuspended in a minimal volume of water. 
P 212 ] For the "ligation" library a DNA (UPPERCASE)-RNA 
(lowercase) hybrid 5-adapter 
( 5 -TACTAATACG ACTCACTaaa- 3 Dharmacon # P- 
002046-01-05) was ligated to the 3-adapted RNA, re- 
verse transcribed with "EcoRI-RT": 
( 5 - G ACTAG CTG G AATTC AAG G ATG CG GTTAAA- 3 ) , PCR 
amplified with two external primers essentially as in El- 
bashir et al 2001 except that primers were "EcoRI-RT" and 



"Pstl 

Fwd " ( 5 - C AG CC AACG CTG C AG ATACG ACTC ACTAAA- 3 ) . 
This PCR product was used as a template for a second 
round of PCR with one hemispecific and one external 
primer or with two hemispecific primers. 
[0213] For the "One tailed" library the 3-Adapted RNA was an- 
nealed to 20pmol primer "EcoRI RT" by heating to 70C and 
cooling O.lC/sec to 30C and then reverse transcribed with 
Superscript II RT (According to instructions, Invitrogen) in 
a 20I volume for 10 alternating 5 minute cycles of 37C 
and 45C. Subsequently, RNA was digested with II 2M 
NaOH, 2mM EDTA at 65C for 10 minutes. cDNA was 
loaded on a polyacrylamide gel, excised and gel-purified 
from excess primer as above (invisible, judged by primer 
run alongside) and resuspended in 131 of water. Purified 
cDNA was then oligo-dC tailed with 400U of recombinant 
terminal transferase (Roche molecular biochemicals), II 
100M dCTP, II 15mM CoCI2, and 41 reaction buffer, to a 
final volume of 201 for 15 minutes at 37C. Reaction was 
stopped with 21 0.2M EDTA and 151 3M NaOAc pH 5.2. 
Volume was adjusted to 1501 with water, Phenol: Bro- 
mochloropropane 10:1 extracted and subsequently pre- 
cipitated with glycogen and 3 volumes of Ethanol. C-tailed 



cDNA was used as a template for PCR with the external 
primers 

"T3-PstBsg(G/l)18"(5 - AATTAACC CTC ACT AAAG G CTG C AG 
GTGCAGGIGGGIIGGGIIGGGIIGN-3 where I stands for Ino- 
sine and N for any of the 4 possible deoxynucleotides), 
and with "EcoRI 

Nested"(5-GGAATTCAAGGATGCGGTTA-3). This PCR 
product was used as a template for a second round of PCR 
with one hemispecific and one external primer or with two 
hemispecific primers. 
[0214] Hemispecific primers were constructed for each predicted 
GAM RNA oligonucleotide by an in-house program de- 
signed to choose about half of the 5 or 3 sequence of the 
GAM RNA corresponding to a TM of about 30-34C con- 
strained by an optimized 3 clamp, appended to the 
cloning adapter sequence (for "One-tailed" libraries 
5-GGNNGGGNNG on the 5 end of the GAM RNA , or TT- 
TAACCGCATC-3 on the 3 end of the GAM RNA. For "Liga- 
tion" libraries the same 3 adapter and 5-CGACTCACTAAA 
on the 5 end). Consequently, a fully complementary 
primer of a TM higher than 60C was created covering only 
one half of the GAM RNA sequence permitting the unbi- 
ased elucidation by sequencing of the other half. 



[0 21 5] CONFIRMATION OF GAM OLIGONUCLEOTIDE SEQUENCE 
AUTHENTICITY OF PCR PRODUCTS: 

[0216] SOUTHERN BLOT: PCR-product sequences were confirmed 
by southern blot (Southern E.M., Biotechnology, 
1992,24:122-39 (1975)) and hybridization with DNA 
oligonucleotide probes synthesized as complimentary 
(antisense) to predicted GAM RNA oligonucleotides. Gels 
were transferred onto a Biodyne PLUS 0.45m, (Pall) posi- 
tively charged nylon membrane and UV cross-linked. Hy- 
bridization was performed overnight with DIG-labeled 
probes at 420C in DIG Easy-Hyb buffer (Roche). Mem- 
branes were washed twice with 2xSSC and 0.1% SDS for 10 
minutes at 420C and then washed twice with 0.5xSSC and 
0.1% SDS for 5 min at 420C. The membrane was then de- 
veloped by using a DIG luminescent detection kit (Roche) 
using anti-DIG and CSPD reaction, according to the man- 
ufacturer's protocol. All probes were prepared according 
to the manufacturers (Roche Molecular Biochemicals) pro- 
tocols: Digoxigenin (DIG) labeled antisense transcripts 
was prepared from purified PCR products using a DIG RNA 
labeling kit with T3 RNA polymerase. DIG labeled PCR was 
prepared by using a DIG PCR labeling kit. 3-DIG-tailed 
oligo ssDNA antisense probes, containing DIG-dUTP and 



dATP at an average tail length of 50 nucleotides were pre- 
pared from lOOpmole oligonucleotides with the DIG 
Oligonucleotide Labeling Kit. 

[O 217 ] CLONE-SEQUENCING: PCR products were inserted into 
pGEM-T (Promega) or pTZ57/T (MBI Fermentas), heat- 
shock transformed into competent JM109 E. coli 
(Promega) and sown on LB-Ampiccilin plates with IPTG 
and Xgal. White and light-blue colonies were transferred 
to duplicate gridded plates, one of which was blotted onto 
a membrane (Biodyne Plus, Pall) for hybridization with DIG 
tailed oligo probes (according to instructions, Roche) 
complimentary to the expected GAM. Plasmid DNA from 
positive colonies was sequenced. 

[0 21 8] LIGATION-PCR DIAGNOSTIC METHOD: To further validate 
predicted GAM PCR product sequence derived from hemi- 
primers, a PCR based diagnostic technique was devised to 
amplify only those products containing also at least two 
additional nucleotides of the non hemi-primer defined 
part of the predicted GAM RNA oligonucleotide. In 
essence, a diagnostic primer was designed so that its 3 
end, which is the specificity determining side, was identi- 
cal to the desired GAM RNA oligonucleotide, 2-10 nu- 
cleotides (typically 4-7, chosen for maximum specificity) 



further into its 3 end than the nucleotide stretch primed 
by the hemi-primer. The hemi-primer PCR product was 
first ligated into a T-cloning vector (pTZ57/T or pGEM-T) 
as described hereinabove. The ligation reaction mixture 
was used as template for the diagnostic PCR under strict 
annealing conditions with the new diagnostic primer in 
conjunction with a general plasmid-homologous primer, 
resulting in a distinct ~200 base-pair product. This PCR 
product can be directly sequenced, permitting the eluci- 
dation of the remaining nucleotides up to the 3 of the 
mature GAM RNA oligonucleotide adjacent to the 3 
adapter. Alternatively, following analysis of the diagnostic 
PCR reaction on an agarose gel, positive ligation reactions 
(containing a band of the expected size) were transformed 
into E. coli. Using this same diagnostic technique and as 
an alternative to screening by Southern-blot colony- 
hybridization, transformed bacterial colonies were 
screened by colony-PCR (Gussow,D. and Clackson,T, Nu- 
cleic Acids Res. 17: 4000 (1989)) prior to plasmid purifi- 
cation and sequencing. 
[0219] Reference is now made to Fig. 13B which is a table sum- 
marizing laboratory results which validate the expression 
of 8 known human miRNA oligonucleotides and 43 novel 



GAM oligonucleotides detected by the BIOINFORMATIC 
OLIGONUCLEOTIDE DETECTION ENGINE 100. The table 
gives additional information on the above GAM oligonu- 
cleotides and contains the following fields: NUMBER: 
refers to the hybridization picture number of Fig. 13A; 
NAME: indicates the known MIR name for the reference 
set, or the GAM's name as given by RosettaGenomics 
nomenclature method; SEQUENCE: 5' to 3' sequence of the 
mature, 'diced' oligonucleotide; SEQUENCED: '+' indicates 
the additional validation of the GAM RNA sequence by se- 
quencing procedure as described hereinabove with refer- 
ence to Fig. 13A. 
[0220] Reference is now made to Fig. 14A, which is a schematic 
representation of a novel human GR polynucleotide herein 
designated GR12731 (RosettaGenomics Ltd. nomencla- 
ture), located on chromosome 9, comprising 2 known hu- 
man MIR oligonucleotides - MIR24 and MIR23, and 2 
novel GAM oligonucleotides, herein designated GAM22 
and GAM116, all marked by solid black boxes. Fig. 14A 
also schematically illustrates 6 non-GAM hairpin se- 
quences, and one non-hairpin sequence, all marked by 
white boxes, and serving as negative controls. By "non- 
GAM hairpin sequences" is meant sequences of a similar 



length to known MIR PRECURSOR sequences, which form 
hairpin secondary folding pattern similar to MIR PRECUR- 
SOR hairpins, and yet which are assessed by the BIOIN- 
FORMATIC OLIGONUCLEOTIDE DETECTION ENGINE 100 
not to be valid GAM PRECURSOR hairpins. It is appreciated 
that Fig. 14A is a simplified schematic representation, re- 
flecting only the order in which the segments of interest 
appear relative to one another, and not a proportional 
distance between the segments. 

[0221] Reference is now made to Fig. 14B, which is a schematic 
representation of secondary folding of each of the MIRs 
and GAMs of GR GR12731 MIR24, MIR23, GAM22 and 
GAM116, and of the negative control non-GAM hairpins, 
herein designated N2, N3, N116, N4, N6 and N7. NO is a 
non-hairpin control, of a similar length to that of known 
MIR PRECURSOR hairpins. It is appreciated that the nega- 
tive controls are situated adjacent to and in between real 
MIR oligonucleotides and GAM predicted oligonucleotides 
and demonstrates similar secondary folding patterns to 
that of known MIRs and GAMs. 

[0222] Reference is now made to Fig. 14C, which is a picture of 
laboratory results of a PCR test upon a YM100 size- 
fractionated "ligation'-library, utilizing a set of specific 



primer pairs located directly inside the boundaries of the 
hairpins. Due to the nature of the library the only PCR am- 
plifiable products can result from RNaselll type enzyme 
cleaved RNA, as expected for legitimate hairpin precursors 
presumed to be produced by DROSHA (Lee et al, Nature 
425 415-419, 2003). Fig 14C demonstrates expression of 
hairpin precursors of known MIR oligonucleotides - 
MIRhsa-23 and MIRhsa-24, and of novel bioinformatically 
detected GAM22 and GAM116 hairpins predicted bioinfor- 
matically by a system constructed and operative in accor- 
dance with a preferred embodiment of the present inven- 
tion. Fig. 14C also shows that none of the 7 controls (6 
hairpins designated N2, N3, N23, N4, N6 and N7 and 1 
non-hairpin sequence designated NO) were expressed. 
N 116 is a negative control sequence partially overlapping 
GAM116. 

[0223] | n the picture, test lanes including template are desig- 
nated "+" and the control lane is designated The con- 
trol reaction contained all the components of the test re- 
action except library template. It is appreciated that for 
each of the tested hairpins, a clear PCR band appears in 
the test ("+") lane, but not in the control ("-") lane. 

[0224] pigs. 14A through 14C, when taken together validate the 



efficacy of the bioinformatic oligonucleotide detection en- 
gine in: (a) detecting known MIR oligonucleotides; (b) de- 
tecting novel GAM PRECURSOR hairpins which are found 
adjacent to these MIR oligonucleotides, and which despite 
exhaustive prior biological efforts and bioinformatic de- 
tection efforts, went undetected; (c) discerning between 
GAM (or MIR) PRECURSOR hairpins, and non-GAM hair- 
pins. 

[0225] it j S appreciated that the ability to discern GAM-hairpins 
from non-GAM-hairpins is very significant in detecting 
GAM oligonucleotides since hairpins are highly abundant 
in the genome. Other MIR prediction programs have not 
been able to address this challenge successfully. 

[0226] Reference is now made to Fig. 15A which is an annotated 
sequence of an EST comprising a novel GAM oligonu- 
cleotides detected by the oligonucleotide detection sys- 
tem of the present invention. Fig. 15A shows the nu- 
cleotide sequence of a known human non-protein coding 
EST (Expressed Sequence Tag), identified as EST72223. 
The EST72223 clone obtained from TIGR database 
(Kirkness and Kerlavage, 1997) was sequenced to yield the 
above 705bp transcript with a polyadenyl tail. It is appre- 
ciated that the sequence of this EST comprises sequences 



of one known miRNA oligonucleotide, identified as hsa- 
MIR98, and of one novel GAM oligonucleotide referred to 
here as GAM25, detected by the BIOINFORMATIC 
OLIGONUCLEOTIDE DETECTION ENGINE 100 (Fig. 2) of the 
present invention. 
[0227] The sequences of the precursors of the known MIR98 and 
of the predicted GAM25 precursors are marked in bold, 
the sequences of the established miRNA 98 and of the 
predicted miRNA-like oligonucleotide GAM25 are under- 
lined. 

[0228] Reference is now made to Figs. 15B, 15C and 15D that are 
pictures of laboratory results, which when taken together 
demonstrate laboratory confirmation of expression of the 
bioinformatically detected novel oligonucleotide of Fig. 
15A. In two parallel experiments, an enzymatically syn- 
thesized capped, EST72223 RNA transcript, was incubated 
with Hela S100 lysate for 0 minutes, 4 hours and 24 
hours. RNA was subsequently harvested, run on a dena- 
turing polyacrylamide gel, and reacted with either a 102nt 
antisense MIR98 probe or a 145nt antisenseGAM25 pre- 
cursor transcript probe respectively. The Northern blot re- 
sults of these experiments demonstrated processing of 
EST72223 RNA by Hela lysate (lanes 2-4, in 15B and 15C), 



into ~80bp and ~22bp segments, which reacted with the 
MIR98 precursor probe (15B), and into ~100bp and ~24bp 
segments, which reacted with the GAM25 precursor probe 
(15C). These results demonstrate the processing of 
EST72223 by Hela lysate into MIR98 precursor and GAM 2 5 
precursor. It is also appreciated from Fig. 15C (lane 1) that 
Hela lysate itself reacted with the GAM25 precursor probe, 
in a number of bands, including a ~100bp band, indicat- 
ing that GAM25-precursor is endogenously expressed in 
Hela cells. The presence of additional bands, higher than 
lOObp in lanes 5-9 probably corresponds to the presence 
of nucleotide sequences in Hela lysate, which contain the 
GAM25 sequence. 
[0229] | n addition, in order to demonstrate the kinetics and 

specificity of the processing of MIR98 and GAM25 precur- 
sors into their respective mature, 'diced' segments, tran- 
scripts of MIR98 and of the bioinformatically predicted 
GAM25 precursors were similarly incubated with Hela 
S100 lysate, for 0 minutes, 30 minutes, 1 hour and 24 
hours, and for 24 hours with the addition of EDTA, added 
to inhibit Dicer activity, following which RNA was har- 
vested, run on a polyacrylamide gel and reacted with 
MIR98 and GAM25 precursor probes. Capped transcripts 



were prepared for in-vitro RNA cleavage assays with 17 
RNA polymerase, including a mZGCS'JpppCS'JG-capping re- 
action using the T7-mMessage mMachine kit (Ambion). 
Purified PCR products were used as template for the reac- 
tion. These were amplified for each assay with specific 
primers containing a 17 promoter at the 5 end and a T3 
RNA polymerase promoter at the 3end. Capped RNA tran- 
scripts were incubated at 30C in supplemented, dialysis 
concentrated, Hela S100 cytoplasmic extract (4C Biotech, 
Seneffe, Belgium). The Hela S100 was supplemented by 
dialysis to a final concentration of 20mM Hepes, lOOmM 
KCI, 2.5mM MgCI2 , 0.5mM DTT, 20% glycerol and pro- 
tease inhibitor cocktail tablets (Complete mini Roche 
Molecular Biochemicals). After addition of all components, 
final concentrations were lOOmM capped target RNA, 
2mM ATP, 0.2mM GTP, 500U/ml RNasin, 25g/ml creatine 
kinase, 25mM creatine phosphate, 2.5mM DTT and 50% 
S100 extract. Proteinase K, used to enhance Dicer activity 
(Zhang et al., EMBO J. 21, 5875-5885 (2002)) was dis- 
solved in 50mM Tris-HCI pH 8, 5mM CaCI2, and 50% 
glycerol, was added to a final concentration of 0.6 mg/ml. 
Cleavage reactions were stopped by the addition of 8 vol- 
umes of proteinase K buffer (200Mm Tris-Hcl, pH 7.5, 



25m M EDTA, 300mM NaCI, and 2% SDS) and incubated at 
65C for 15min at different time points (0, 0.5, 1, 4, 24h) 
and subjected to phenol/chloroform extraction. Pellets 
were dissolved in water and kept frozen. Samples were 
analyzed on a segmented half 6%, half 13% polyacrylamide 
1XTBE-7M Urea gel. 

[0230] The Northern blot results of these experiments demon- 
strated an accumulation of a~22bp segment which re- 
acted with the MIR98 precursor probe, and of a ~24bp 
segment which reacted with the GAM25 precursor probe, 
over time (lanes 5-8). Absence of these segments when 
incubated with EDTA (lane 9), which is known to inhibit 
Dicer enzyme (Zhang et al., 2002), supports the notion 
that the processing of MIR98 and CAM25 precursors into 
their ' diced' segments is mediated by Dicer enzyme, 
found in Hela lysate. Other RNases do not utilize divalent 
cations and are thus not inhibited by EDTA. The molecular 
sizes of EST72223, MIR-98 and GAM25 and their corre- 
sponding precursors are indicated by arrows. 

[0231] pig. 15D present Northern blot results of same above ex- 
periments with GAM25 probe (24nt). The results clearly 
demonstrated the accumulation of mature GAM25 
oligonucleotide after 24 h. 



[0232] jo validate the identity of the band shown by the lower 
arrow in figs. 15C and 15D, a RNA band parallel to a 
marker of 24 base was excised from the gel and cloned as 
in Elbashir et al (2001) and sequenced. 90 clones corre- 
sponded to the sequence of mature GAM25 oligonu- 
cleotide, three corresponded to GAM25* (the opposite 
arm of the hairpin with a 1-3 nucleotide 3 overhang) and 
two to the hairpin-loop. 

[0233] GAM25 was also validated endogenously by sequencing 
from both sides from a HeLa YM100 total-RNA "ligation" 
libraries, utilizing hemispecific primers as described in 
Fig. 13. 

[0234] Taken together, these results validate the presence and 
processing of a novel MIR-like oligonucleotide, CAM25, 
which was predicted bioinformatically. The processing of 
this novel GAM oligonucleotide product, by Hela lysate 
from EST72223, through its precursor, to its final form 
was similar to that observed for known miRNA oligonu- 
cleotide, MIR98. 

[0235] Transcript products were 705nt (EST72223), 102nt (MIR98 
precursor), 125nt (GAM25 precursor) long. EST72223 was 
PCR amplified with T7-EST 72223 forward primer: 
5 -TAATAC G ACTC ACTATAG G CCCTTATTAG AG G ATTCTG CT 



-3 and T3-EST72223 reverse 

primer: " - AATTAACCCTC ACTAAAG GTT I I I I I I TCCTGAGA 
CAGAGT-3.MIR98 was PCR amplified using EST72223 as a 
template with T7MIR98 forward primer: 
5 -TAATACG ACTCACTATAGGGTG AG GTAGTAAGTTGTATTG 
TT-3and T3MIR98 reverse primer: 

5 -AATTAACCCTCACTAAAGGGAAAGTAGTAAGTTGTATAGT 
T-3.GAM25 was PCR amplified using EST72223 as a tem- 
plate with GAM25 forward primer: 
5 - GAG G C AG G AG AATTG CTTG A- 3 and T3-EST72223 re- 
verse 

p r i m e r : 5 - AATTAACCCTC ACTAAAG G CCTG AG AC AG AGTCT 
TGCTC-3. 

[0236] it is appreciated that the data presented in Figs. 15A, 15B, 
15C and 15D when taken together validate the function of 
the bioinformatic oligonucleotide detection engine 100 of 
Fig. 2. Fig. 15A shows a novel GAM oligonucleotide bioin- 
formatically detected by the BIOINFORMATIC OLIGONU- 
CLEOTIDE DETECTION ENGINE 100, and Figs. 15C and 15D 
show laboratory confirmation of the expression of this 
novel oligonucleotide. This is in accord with the engine 
training and validation methodology described herein- 
above with reference to Fig. 3. 



[0237] Reference is now made to Fig. 16, which presents pictures 
of laboratory results, which demonstrate laboratory con- 
firmation of excision (dicing) of the bioinformatically de- 
tected novel GAM HIV1 miRNA-like oligonucleotides, 
herein designated GAM2032-B, GAM3249-A, GAM 507-B 
and GAM1016-B from their predicted precursors by their 
incubation in HeLa S-100 lysate as described in Fig. 15. 

[0238] Fig.l6A presents the entire 5-UTR of HIV1 (U5R) contain- 
ing two predicted GAM precursor genes, in bold; GAM 
2032 and GAM3249. The bioinformatically predicted ma- 
ture GAMs are depicted with underscore, the 5-most is 
GAM 2032-B and the second is GAM 3249-A. GAM 
2032-B matches the known HIV1 RNA structure named 
TAR to which the TAT protein binds (Nature 1987. 
330:489-93). 

[0239] pig. 16B and Fig. 16C depict Northern blot analysis of 
GAMs in U5R, hybridized with predicted mature GAM 
oligonucleotide probes GAM 2032-B, and GAM 3249-A, 
respectivley. The molecular size of the entire U5R tran- 
script, 355nt, is indicated by arrow. The predicted molec- 
ular sizes of GAM 2032-B, and GAM 3249-A are 22nt and 
17nt respectively. The 22nt molecular marker is indicated 
by arrow. Lanes: 1-Hela lysate; 2-U5R transcript in HeLa 



Lysate without incubation. 3- U5R transcript incubated 
overnight with Hela lysate. 

[0240] pigs. 16D and 16E present partial transcripts of HIV1 RNA 
reacted with predicted mature HIV1-GAM oligonucleotide 
probes. In each figure, the experimental transcript se- 
quence is shown, annotated in bold is the predicted GAM 
precursor, and in underscore the predicted mature GAM. 
Northern blot analysis of GAM precursors for GAM507-B 
(Fig. 16D), and GAM1016-B (Fig. 16E). The transcript sizes 
are 163nt for GAM507-B transcript and 200nt for 
GAM1016-B transcript. The predicted molecular sizes of 
GAM507-B and GAM1016-B are both 24nt. The 22nt 
molecular marker is indicated by arrow. Lanes: 1- Tran- 
script in HeLa Lysate without incubation. 2- Transcript in- 
cubated overnight with HeLa lysate. 

[° 241 ] It is appreciated that the sequence of the reacting bands 
in the foreseen sizes comprise sequences of novel GAM 
genes, referred to here as GAM 2032-B, GAM 3249-B, 
GAM 507-B , GAM 1016-B , detected by the bioinformatic 
gene detection engine 100 of the present invention, de- 
scribed hereinabove with reference to Fig. 2. 
DETAILED DESCRIPTION OF LARGE TABLES 

[0242] Table 1 comprises data relating the SEQ ID NO of oligonu- 



cleotides of the present invention to their corresponding 
GAM NAME, and contains the following fields: GAM SEQ- 
ID: GAM SEQ ID NO, as in the Sequence Listing; GAM 
NAME: Rosetta Genomics Ltd. nomenclature (see below); 
GAM RNA SEQUENCE: Sequence (5 x to 3 % ) of the mature, 
x diced x GAM RNA; GAM ORGANISM: identity of the organ- 
ism encoding the GAM oligonucleotide; GAM POS: Dicer 
cut location (see below); and 

[0243] Table 2 comprises detailed textual description according 
to the description of Fig.l of each of a plurality of novel 
GAM oligonucleotides of the present invention, and con- 
tains the following fields: GAM NAME: Rosetta Genomics 
Ltd. nomenclature (see below); GAM ORGANISM: identity 
of the organism encoding the GAM oligonucleotide; PRE- 
CUR SEQ-ID:GAM precursor Seq-ID, as in the Sequence 
Listing; PRECURSOR SEQUENCE: Sequence (5 X to 3^) of the 
GAM precursor ; GAM DESCRIPTION: Detailed description 
of GAM oligonucleotide with reference to Fig.l; and 

[0244] Table 3 comprises data relating to the source and location 
of novel GAM oligonucleotides of the present invention, 
and contains the following fields: GAM NAME: Rosetta Ge- 
nomics Ltd. nomenclature (see below); PRECUR SEQ-ID: 
GAM precursor SEQ ID NO, as in the Sequence Listing; 



GAM ORGANISM: identity of the organism encodes the 
GAM oligonucleotide; SOURCE: For human GAM- 
chromosome encoding the human GAM oligonucleotide, 
otherwise- accession ID (GenBank,NCBI); STRAND: Orien- 
tation of the strand, '+' for the plus strand, '-' for the mi- 
nus strand; SRC-START OFFSET: Start offset of GAM pre- 
cursor sequence relative to the SOURCE; SRC-END OFFSET: 
End offset of GAM precursor sequence relative to the 
SOURCE; and 

[0245] Table 4 comprises data relating to GAM precursors of 

novel GAM oligonucleotides of the present invention, and 
contains the following fields: GAM NAME: Rosetta Ge- 
nomics Ltd. nomenclature (see below); PRECUR SEQ-ID: 
GAM precursor Seq-ID, as in the Sequence Listing; GAM 
ORGANISM: identity of the organism encoding the GAM 
oligonucleotide; PRECURSOR-SEQUENCE: GAM precursor 
nucleotide sequence (5^ to 3^); GAM FOLDED PRECURSOR 
RNA: Schematic representation of the GAM folded precur- 
sor, beginning 5^ end (beginning of upper row) to 3 X end 
(beginning of lower row), where the hairpin loop is posi- 
tioned at the right part of the draw; and 

[0246] Table 5 comprises data relating to GAM oligonucleotides 
of the present invention, and contains the following fields: 



GAM NAME: Rosetta Genomics Ltd. nomenclature (see be- 
low); GAM ORGANISM: identity of the organism encoding 
the GAM oligonucleotide; GAM RNA SEQUENCE: Sequence 
(5^ to 3^) of the mature, x diced^ GAM RNA ; PRECUR 
SEQ-ID : GAM precursor Seq-ID, as in the Sequence List- 
ing; GAM POS: Dicer cut location (see below); and 

[0247] Table 6 comprises data relating SEQ ID NO of the GAM 

target gene binding site sequence to TARGET gene name 
and target binding site sequence, and contains the follow- 
ing fields: TARGET BINDING SITE SEQ-ID: Target binding 
site SEQ ID NO, as in the Sequence Listing; TARGET OR- 
GANISM: identity of organism encode the TARGET gene; 
TARGET: GAM target gene name; TARGET BINDING SITE 
SEQUENCE: Nucleotide sequence (5^ to 3^) of the target 
binding site; and 

[0248] Table 7 comprises data relating to target-genes and bind- 
ing sites of GAM oligonucleotides of the present inven- 
tion, and contains the following fields: GAM NAME: 
Rosetta Genomics Ltd. nomenclature (see below); GAM 
ORGANISM: identity of the organism encoding the GAM 
oligonucleotide; GAM RNA SEQUENCE: Sequence (5 X to 3 X ) 
of the mature, ^diced^ GAM RNA; TARGET: GAM target 
gene name; TARGET REF-ID: Target accession number 



(GenBank); TARGET ORGANISM: identity of organism en- 
code the TARGET gene; UTR: Untranslated region of bind- 
ing site/s (3" or 5"); TARGET BS-SEQ: Nucleotide sequence 
(5^ to 3^) of the target binding site; BINDING-SITE-DRAW: 
Schematic representation of the binding site, upper row 
represent 5 X to 3 X sequence of the GAM, Lower row rep- 
resent 3 X to 5 X Sequence of the target; GAM POS: Dicer 
cut location (see below); and 
[0249] Table 8 comprises data relating to functions and utilities 
of novel GAM oligonucleotides of the present invention, 
and contains the following fields: GAM NAME: Rosetta Ge- 
nomics Ltd. nomenclature (see below); GAM RNA SE- 
QUENCE: Sequence (5 X to 3 N ) of the mature, x diced x GAM 
RNA; GAM ORGANISM: identity of the organism encoding 
the GAM oligonucleotide; TARGET:GAM target gene name; 
TARGET ORGANISM: identity of organism encode the TAR- 
GET gene; GAM FUNCTION: Description of the GAM func- 
tions and utilities; GAM POS: Dicer cut location (see be- 
low); and 

[0250] Table 9 comprises data of gene function references- Bibli- 
ography and contains the following fields: GAM NAME: 
Rosetta Genomics Ltd. nomenclature (see below); GAM 
RNA SEQUENCE: Sequence (5^ to 3 X ) of the mature, 



^dicecT GAM RNA gene; GAM ORGANISM: identity of the 
organism encoding the GAM oligonucleotide; TARGET: 
GAM target gene name; TARGET ORGANISM: identity of 
organism encode the TARGET gene; REFERENCES: list of 
references relating to the target gene; GAM POS: Dicer cut 
location (see below); and 

[0251] Table 10 comprises data relating to novel GR (Genomic 
Record) polynucleotides of the present invention, and 
contains the following fields: GR NAME: Rosetta Genomics 
Ltd. nomenclature (see below); GR ORGANISM: identity of 
the organism encoding the GR polynucleotide; GR DE- 
SCRIPTION: Detailed description of a GR gene cluster, with 
reference to Fig. 9; and 

[0252] The following conventions and abbreviations are used in 
the tables: The nucleotide 'IT is represented as T in the 
tables, and 

[0253] gam NAME or GR NAME are names for nucleotide se- 
quences of the present invention given by RosettaGe- 
nomics Ltd. nomenclature method. All GAMs/GRs are des- 
ignated by GAMx/GRx where x is a unique ID. 

[0254] GAM POS is a position of the GAM RNA on the GAM PRE- 
CURSOR RNA sequence. This position is the Dicer cut lo- 
cation: A indicates a probable Dicer cut location; B indi- 



cates an alternative Dicer cut location. 
[0255] All human nucleotide sequences of the present invention 
as well as their chromosomal location and strand orienta- 
tion are derived from sequence records of UCSC-hgl6 
version, which is based on NCBI, Build34 database (April, 
2003). 

[0256] All viral sequences of the present invention as well as their 
genomic location are derived from NCBI, RefSeq database. 



